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ABSTRACT 


Self-consistent  computer  calculations  with  the 
McMillan  tunneling  model  of  the  proximity  sandwich  have 
been  carried  out  and  a  detailed  analysis  of  the  tunneling 
density  of  states  has  been  made.  Also  tunneling  experiments 
have  been  performed  into  the  aluminum  side  of  some  aluminum- 
tin  proximity  sandwiches  and  into  the  tin  side  of  others 
evaporated  at  room  temperature  onto  an  oxidized  aluminum 
electrode  B.  The  coupling  between  the  aluminum  and  tin 
proximity  films  was  made  weak  by  purposely  allowing  slight 
oxidation  to  occur  at  the  interface.  The  properties  of  the 
sandwich  depend  on  the  amount  of  this  oxidation.  When  B  is 
superconducting,  the  double-peaked  structure  of  the  McMillan 
density  of  states  is  resolved  in  the  tunneling  conductance 
as  a  function  of  voltage  and  the  range  of  its  observability 
has  been  established.  When  B  is  normal,  the  normalized 
tunneling  conductance  between  B  and  the  aluminum  side  of  the 
sandwich  in  the  vicinity  of  the  critical  temperature  of  the 
sandwich  is  markedly  different  from  that  of  junctions  formed 
between  B  and  a  BCS  superconductor.  The  normalized 
conductance  between  B  when  it  is  normal  and  the  tin  side 
of  the  sandwich  in  the  vicinity  of  the  critical  temperature 
of  the  sandwich  is  about  the  same  as  that  of  a  junction 
formed  between  B  and  a  BCS  superconductor.  The  comparison 


■ 


of  these  experimental  results  shows  satisfactory 
quantitative  agreement  with  the  predictions  of  the  McMillan 
model . 
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CHAPTER  I 


INTRODUCTION 


Within  the  past  few  years  there  has  been  considerable 
interest  in  experiments  on  the  proximity  effect  between 
normal  (N)  and  superconducting  (S)  materials,  usually 
prepared  by  evaporation. 

Theoretical  difficulties  arise  in  calculating  the 
superconducting  properties  of  proximity  sandwiches,  because 
in  such  systems  the  order  parameter  is  space  dependent. 
Solutions  have  been  obtained  only  in  the  vicinity  of  the 
critical  temperature,  T  .  The  De  Gennes  -  Werthamer  theory, 
which  is  representative  of  this  approach,  has  been  quite 
well  verified  in  the  range  of  its  validity  by  a  number  of 
experiments . 

McMillan  avoids  the  problem  of  a  space  dependent 
order  parameter  by  assuming  that  the  N  and  S  films  are 
weakly  coupled  and  that  each  metal  is  sufficiently  thin 
that  the  superconducting  properties  are  constant  across 
it.  He  is  then  able  to  obtain  a  complete  solution  for  his 
model  valid  for  all  temperatures  T  <  T^ .  A  few  experimenters 
found  qualitative  agreement  with  the  McMillan  model,  but 
their  experiments  did  not  quite  fulfil  the  assumptions 


of  the  model. 
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The  purpose  of  this  theses  is  to  test  the  validity 
of  the  McMillan  model  for  aluminum-tin  proximity  sandwiches. 

We  have  carried  out  self  -  consistent  computer 
calculations  of  the  McMillan  model  to  obtain  a  detailed 
analysis  of  the  tunneling  density  of  states.  The  results 
for  Al-Sn  sandwiches  are  presented  in  Chapter  III  in 

graphical  form. 

The  tunneling  experiments  have  been  performed  into 
both  sides  of  Al-Sn  proximity  sandwiches  evaporated  at 
room  temperature  onto  an  oxidized  aluminum  electrode  which 
may  be  either  normal  or  superconducting.  The  coupling 
between  the  N  and  S  films  was  made  weak  by  allowing  slight 
oxidation  to  occur  at  the  interface.  Comparison  of  the 
experimental  tunneling  conductivity  with  the  calculated  one 
shows  satisfactory  quantitative  agreement  and  the  results 
are  presented  in  Chapter  VI. 

Chapter  II  contains  a  brief  review  of  previous 
work  done  on  the  proximity  effect.  In  Chapters  IV  and  V 
the  technology  of  the  tunnel  junctions  preparation,  the 
apparatus  used  and  the  mean  free  path  are  discussed. 

The  results  of  Chapter  III. 2  and  VI. 2. a,  b,  and  c 
will  be  published  in  the  Physical  Review  (Vrba  J. ,  Woods 
S . B. ,  April , 1971) . 

Unless  otherwise  stated, the  Gaussian  units  are  used 


in  the  theses. 
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CHAPTER  II 

THE  REVIEW  OF  PREVIOUS  WORK  ON  PROXIMITY  EFFECT 

II. 1.  General  Considerations 

When  a  superconductor  is  placed  in  very  close  contact 
with  a  normal  metal,  the  properties  of  both  are  changed.  The 
normal  metal  may  become  superconducting,  while  the 
transition  temperature  and  gap  of  the  superconductor  are 
lowered.  If  two  different  superconductors  are  placed  into 
proximity  contact,  their  properties  change  and  the  resulting 
gap  and  critical  temperature  are  somewhere  in  between  the 
values  of  each  superconductor. 

The  early  experimental  works  of  Meissner  (1960)  on 
contacts  between  crossed  wires  of  tin  plated  with  different 
metals,  Rose-Innes  and  Serin  (1961)  on  gold  plates  coated 
by  evaporated  tin.  Smith  et  al  (1961)  on  evaporated  films 
of  lead  and  silver,  and  of  Simmons  and  Douglass  (1962)  on 
films  of  tin  and  silver  aroused  the  interest  of  other 
investigators  in  the  proximity  effect. 

We  may  divide  the  attempts  to  solve  the  problem 
of  the  proximity  effect  into  two  main  groups: 

(a)  works  dealing  with  the  critical  temperature  of  the 
layered  structure  only,  and 

(b)  works  dealing  with  the  excitation  spectrum  of  the 
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layered  structure. 

The  smallness  of  the  order  parameter  A  around  the  critical 

temperature  T  has  been  employed  in  theoretical 

c 

considerations  included  in  group  (a)  .Therefore  the  results 

from  group (a ) cannot  be  used  at  lower  temperatures.  Group  (b) 

includes  both  works  valid  around  T  only  and  works  valid 

c 

at  all  temperatures  T  <  T^ . 

Common  to  the  theoretical  approach  to  the  proximity 
effect  is  the  notion  that  superconductivity  in  bulk  material 

is  characterised  by  an  electron  pair  correlation  over  a 

-4  -5 

distance  of  the  order  of  10  to  10  cm  and  that  this 
correlation  should  extend  a  similar  distance  into  a  normal 
metal  in  contact  with  a  superconductor.  Deutscher  and 
De  Gennes  (1969)  show  that  if  normal  metal  is  "clean",  i.e., 
the  mean  free  path  is  large  compared  to  the  coherence 
length  £  ,  then  the  probability  amplitude  F  =  <  >  for 

finding  a  Cooper  pair  at  a  distance  |x|  from  the  N-S 
boundary  has  the  form 

F=$(X)e  for  |x|  large,  II. 1 

where  0 (X)  is  a  slowly  varying  function  of  X  and 
K*'1  =  -ftv  2-rrk^T,  v^.,  is  the  Fermi  velocity  in  the  N-film 
and  T  is  the  temperature.  When  T  goes  to  zero  then,  as 
pointed  out  by  Falk  (1963)  ,  the  decrease  of  F  with  X  is 
very  slow,  approximately  as  1/X.  If  the  normal  metal  is 
dirty  (  c  )  ,  the  leakage  of  Cooper  pairs  is  ruled  by 

a  diffusion  process  and  the  exponential  form  II. 1  still 

-1  1/2 
holds,  but  K  is  replaced  by  £N  =  ("hvFN  ^  /  6iTkBT) 
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II. 2.  Critical  Temperature  of  Proximity  Sandwich. 

A  simple  picture  of  the  proximity  effect  was  presented 
by  Cooper  (1961,  1962).  He  considers  a  normal  metal  of 
thickness  d^  in  perfect  contact  (  no  oxide  barrier  )  with  a 
superconductor  of  thickness  dg,  both  with  the  same  Fermi 
energy  and  the  same  electronic  effective  mass.  Due  to  the 
coherence  distance  associated  with  zero  momentum  pairs,  the 
superconducting  correlation  can  extend  into  the  volume  of 
normal  metal,  so  that  a  pair  samples  an  effective  attractive 
potential  which  is  simply  the  spatial  average  of  the 
potential  on  the  two  sides  of  the  interface.  Assuming  that 
both  d^T  and  d0  are  smaller  than  the  coherence  lengths  and 
that  [N(0)V]n=0  in  the  normal  metal  and  [N(0)V]s^0  in  the 
superconductor,  Cooper  writes  for  the  effective  BCS  coupling 
constant  of  the  proximity  sandwich 


[N(0)V]eff 


[N  ( 0 )  V]  s  . 


II. 2 


The  BCS  relationship  for  the  resulting  critical  temperature 
of  the  sandwich  may  be  written 


1. 14H(jj  e 


[N(°)V]eff  _ 


However  in  actual  experiments  some  adsorbed  gas  is  always 
present  in  between  the  two  metals.  These  atoms  form 
a  barrier  which  tends  to  separate  metals.  Therefore  one 
should  expect  that  eq.II.2  gives  the  greatest  possible 


6 


reduction  of  the  effective  electron  -  electron  interaction. 
Cooper  suggests  that  the  actual  reduction  factor  should  have 
the  form 


+  3d 


N 


0  = 


6  5 


t 


where  3  is  determined  by  the  barrier  between  the  films. 

A  more  detailed  theory  was  developed  by  De  Gennes  and 
Guyon  (1963)  and  De  Gennes  (1964a).  They  restricted  their 
considerations  to  dirty  films  in  which  the  mean  free  path  £ 
is  small  compared  to  the  coherence  length  6.  They  expect 
to  get  better  agreement  with  theory  for  dirty  systems 
(i  <  6)  than  for  clean  ones,  because  the  detailed  atomic 
situation  at  the  interface  in  a  dirty  system  is  not 
important  and  large  scale  motion  of  superconducting 
electrons  is  ruled  by  the  diffusion  equation. 

De  Gennes  (1964a)  assumes  that  the  electrons  are 
coupled  by  a  point  interaction,  V (r^) 6  (r^-r ^ ) .  In  the 
S-film  he  puts  V=Vg  positive  and  in  the  N-film  V=VN 
may  have  either  sign.  The  interaction  is  cut  off  at 
frequency  or  to^.  He  restricts  all  considerations  to 
the  vicinity  of  the  transition  temperature  T  of  the 
proximity  sandwich  and  assumes  that  the  pair  potential 
A(r)  is  small  at  all  points  in  space  when  the  temperature 
is  close  to  the  transition  point  (  this  arises  from  the 
assumption  that  the  superconductivity  transition  is  of  the 
second  order).  Then  he  treats  A(r)  as  a  perturbation  and 


is  able  to  write  a  linear  integral  equation  for  A (in  units 
'fi  =  k  =  c  =  1)  : 
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A  (r)  =  V  (r)  £ 

03 

H  (r , r ' )  =T  £ 

03  ~  ~ 

03 


/  d3r'A(r')Hu(r,r-)  , 


V  1  1  x  + 

£  ■ — — ■; —  — : —  w  (r)w 

c  -103  e  +103  n  ~  m 
n  ,m  n  m 


(r )  w 
~  n 


(r  '  )  w  (r '  ) 
m  ~ 


II. 3 

II.  4 


where  03=2ttT  (v+1/2),  £  represents  a  sum  over  all  (positive 

or  negative  or  zero)  integers  v  and  T  is  the  temperature. 

The  functions  w  are  one-electron  wave 

n 

functions  in  the  normal  state  defined  by 


El  + 

2m 


U (r )  ]  w 


n 


e  w 
n  n 


They  include  the  effect  of  impurity  and  boundary  scattering. 
The  boundary  conditions  at  the  interface  are 

A/NV  continuous 


II.  5 


(D/V)dA/dX  continuous 

where  D  =  v^£  and  N  is  the  density  of  states  at  the 

Fermi  energy  in  the  normal  state. 

Applying  the  Cooper  assumption  De  Gennes  is  able  to 
derive  the  effective  BCS  coupling  constant 


[  NV  ]  ff= 


V  M  2  d ,  +  \7  T]3  d 
N  N  N  S  S  S 

NMdM+N cd 
N  N  S  S 


II. 6 


Eq.  II. 6  differs  from  Cooper's  eq .  II. 2  by  the  weighting 

factors  which  are ,  according  to  De  Gennes,  dlTItT  and  d^N„, 

N  N  S  S 

while  according  to  Cooper  they  are  d^  and  dg . 
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The  early  experiments  by  Smith  et  al  (1961)  on  Pb  -  Ag 
systems  and  by  Simmons  and  Douglass  (1962)  on  Sn  -  Ag  systems 
might  be  interpreted  in  terms  of  Cooper's  approximation, 
but  the  mean  free  path  and  coherence  length  were  not 
controlled  and  therefore  the  assumptions  for  approximation 
cannot  be  verified. 

De  Gennes  makes  the  following  approximation  for  N-S 
sandwiches  with  thicknesses  d  and  d  larger  than  the 

o  IM 

corresponding  coherence  lengths  and  6^:  in  the  kernel 

H  (r,r')  (eq.  II. 4)  he  retains  only  the  lowest  frequency 

component  H  (coo  =  ±ttt)  for  which  the  kernel  has  maximum 
oo  o 

range.  All  other  components  are  replaced  by  suitable 
normalized  6  functions.  He  is  then  able  to  derive  the 
critical  temperature  of  the  sandwich  (in  units  -fi=k=c=l)  , 


TC  =  Tlira  +  "DS- 


~  -2Kd„ 

-  e  N  , 

3 


II. 7 


(d0+b  ) 
v  S  00 


where  critical  temperature  when  d^-*-00 ,  extrapolation 

length  b  =  DgN^coth (Kd^) /KD^N^ ,  b^  =  lim  b  and 

dN^°° 


2  2ttT  ,  ,  2NNVN  .  ^  TX  _  ..  .  , 

K  =  -rr-  (  l"T~7^rT7~  )•  Ec3*  H*7  was  first  proposed  on  a 


D 


N 


purely  experimental  basis  by  Hilsch  (1962)  for  Cu  -  Pb 
sandwiches . 

If  there  was  a  thin  oxide  layer  between  the  N  and  S 
films,  the  temperature  shift  Tc  -  T^m  would  be  reduced 
but  the  value  of  K  would  remain  the  same  (De  Gennes) . 


■ 
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Thus  K  is  a  characteristic  length  which  can  be  measured 
even  in  the  presence  of  impurities  in  the  interface  (  as 
long  as  the  impurity  content  is  somehow  kept  unchanged 
from  the  sample  to  sample  ) . 

Werthamer  (1963)  is  able,  by  transforming  the  integral 
equation  II. 3  into  a  differential  equation  and  assuming 
that  the  N  -  side  has  a  small  (negligible)  critical 
temperature  TCN ,  to  derive  equations  for  the  critical 
temperature  of  a  superposed  sandwich.  The  refined 
version  by  Hauser,  Theuerer  and  Werthamer  (1964)  leads  to 


NS^kStan(kSdS)  =  VNkNtanh(kNdN}' 


x(?sks)  =  In 


CS 


“X("CNkN}  =  ln 


T 


II.  8 


CN 


arir?  v  ^  7  \  =  W  M  /9  4-7  -  W  H  /?  ) 

A.  '  S-J  '  *  N  /  —  —/  “  /  *  \  -*-/  / 


where  V  is  digamma  function  (see  e.g.  Davis  (1935)),  T 

O 

is  the  critical  temperature  of  the  S  -  film  of  the  sandwich, 
and  K  =  -ftv  £ /6tt kDT„  plays  the  role  of  coherence  length. 

Hauser  et  al  (1964)  used  eqs.  II. 8  to  interpret  their 
data  on  the  critical  temperature  of  Cu  -  Pb  and  Pt  -  Pb 
sandwiches,  and  Hauser  and  Theuerer  (1965)  tested  the 
validity  of  eqs.  II. 8  for  two  superconductors  in  proximity 
(  Pb  -  Al  sandwich  ) .  In  both  cases  the  agreement  with 
theory  was  satisfactory.  Hauser  et  al  (1966)  also  extended 
eqs.  II. 8  for  the  superposed  films  of  superconductor  and 
magnetic  film  (  using  the  Abrikosov  -  Gor ' kov  model  of 
superconductivity  in  dilute  magnetic  alloys  )  and  they 
found  a  reasonable  agreement  with  experiment. 
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I I . 3  Excitation  Spectrum  of  Proximity  Sandwich. 

Fulde  and  Maki  (1965) ,  in  order  to  explain  the 
experiments  of  Woolf  and  Reif  (1965)  on  superimposed  films 
of  superconductor  and  paramagnetic  material,  calculated  the 
density  of  states  in  superimposed  sandwiches  of  both 
superconductor  and  normal  nonmagnetic  metal  and  superconductor 
and  magnetic  metal.  They  obtained  the  expression  for  the 
density  of  states  N(r,w)  valid  from  both  sides  of  the 
superimposed  sandwich, 

|  A (r)  |  2  2  2 

N(r,m)  =  N  ( 0 )  [1+  — -  JO-Z® -  ]  ,  ii. 9 

2  (m2+a2)2 


where  ct  =  x  k2/6,  is  the  transport  collision  time, 

and  k  =  k  (or  k  )  is  determined  by  eqs.  II. 8.  In  case  of 

o  IN 

a  contact  between  a  superconducting  and  a  paramagnetic 
film  the  spin  -  flip  scattering  relaxation  constant  due 
to  the  paramagnetic  ions  is  introduced  and  k^  in  eq.II.8 
is  replaced  by  k^  ,  for  which  new  equation  is  valid  : 


-x( 


ITT  T_ 

s  C 


-  ^kp>  =  ln  T 


c 

Cp 


Expression  II. 9  is  equivalent  to  the  one  obtained  by  De 
Gennes  (1964b)  in  the  gapless  region  for  superconductors 
with  paramagnetic  impurities  or  type  II  superconductors 
with  a  characterising  the  strength  of  the  depairing  effect. 
De  Gennes  and  Mauro  (1965)  pointed  out  that  Fulde  and 


Maki ' s  derivation  of  eq .  II. 9  is  valid  only  for  homogeneous 
samples.  They  showed  that  the  expansion  of  the  density-of- 
states  in  powers  of  the  order  parameter  A(rJ  does  not 
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converge,  but  a  similar  expansion  for  the  tunneling 
characteristic  I(V)  does  converge  at  sufficiently  high 
temperatures.  Then  Fulde  and  Maki's  approximation  is  valid 


(K-v^r1 


film  and  for  dAT  > 

N 


for  the  N-film. 

Their  results  were  generalized  and  tested  experimentally 
by  Guyon  et  al  (1966)  for  superimposed  sandwiches  of  dirty 
films  of  InBi  alloys  and  Zn.  The  agreement  with  the  theory 
is  very  good.  Claeson  and  Gygax  (1966)  and  Claeson  et  al 
(1967)  investigating  the  tunneling  into  Pb-Al,  Pb-Ag ,  and 
Pb-Mn  sandwiches  and  Hideki  Tsuya  (1967)  investigating 
Ag-Sn,  Ag-Pb,  Fe-Sn,  and  Fe-Pb  sandwiches  found  general 
agreement  with  the  Fulde  and  Maki  calculations. 

McMillan  (1968)  circumvented  the  problems  with  the 
solution  of  integral  equation  II. 3  for  the  pair  potential  A 
by  proposing  a  simple,  soluble  tunneling  model  of  proximity 
effect.  The  potential  barrier  is  assumed  to  separate  the 
N  and  S  films  and  tunneling  through  this  barrier  is  treated 
by  the  tunneling  Hamiltonian  method.  Using  this  model 
McMillan  was  able  to  calculate  the  tunneling  density  of 
states  in  each  film  of  an  N  -  S  sandwich  for  all  temperatures 
T  <  T  .  The  thickness  of  each  film  is  assumed  to  be  smaller 
than  the  corresponding  coherence  length  and  therefore  the 
calculation  does  not  involve  the  coherence  length  directly. 
Kaiser  and  Zuckermann  (1970)  extended  the  McMillan  model 
to  the  case  when  the  normal  metal  contains  magnetic 
impurities.  They  also  made  a  slight  numerical  correction 


. 
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to  one  relation  among  the  model  parameters. 

Freake  and  Adkins  (1969)  measured  the  tunneling 
conductance  between  normal  metal  and  Cu-Pb  sandwiches  at 
temperatures  of  0.06°  -  0.08°K.  Adkins  and  Kington  (1969) 
measured  tunneling  from  normal  aluminum  into  Ag-Pb, 

Cu-Pb,  and  Sn-Pb  sandwiches . Both  groups  of  authors  found 
qualitative  agreement  with  McMillan's  predictions,  which  is 
as  much  as  could  be  expected  since  their  films  were  rather 
strongly  coupled,  and  the  theory  only  applies  strictly  to 
weak  coupling  between  the  two  films. 

II . 4  Experimental  Difficulties. 

To  conclude  this  chapter  the  experimental  difficulties 
connected  with  proximity  effect  experiments  are  discussed. 
All  interpretations  of  the  experimental  results  should  be 
done  with  regard  to  these  problems. 

(a)  As  pointed  out  first  by  Rose-Innes  and  Serin  (1961) , 
the  change  of  critical  temperature  of  a  superconductor 
evaporated  on  top  of  a  normal  metal  may  be  caused  by 
spurious  interdiffusion  effects.  This  effect  was 
studied  in  more  detail  by  Chiou  and  Klokholm  (1964) 
on  Ag-Sn  sandwiches,  Klokholm  and  Chiou  (1964)  on  Au -Sn 
sandwiches,  Caswell  (1964)  on  Pb-Cu,  Pb-Ag ,  and  Pb-Au 
sandwiches  and  by  Van  Gurp  (1963)  on  Sn-Al,  In-Al ,  and 
In-Mg.  It  was  concluded  that  when  intermetallic 
compounds  are  known  to  exist  the  change  of  critical 
temperature  is  determined  by  the  formation  of  compounds 
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rather  than  by  true  proximity  effect. 

(b)  Hilsch  (1962)  first  observed  irreversibility,  i.e.  the 
dependence  of  transition  temperature  of  a  sandwich  on 
the  order  of  deposition  of  the  metals.  Hauser  et  al 
(1966)  argue  on  the  basis  of  their  improved  experimental 
arrangement  that  the  irreversibility  is  caused  by  strong 
dependence  of  resistivity  (or  mean  free  path)  of  the 
first  few  layers  of  a  deposited  film  on  the  type  of 
surface  where  growth  is  initiated.  They  rule  out  the 
possibility  of  a  thin  oxide  layer  at  the  film  interface 
(Hauser  et  al  (1964) ) .  However  on  the  basis  of  our 
experience  gained  during  the  present  work  we  believe, 
that  at  least  for  metals  which  oxidize  readily,  the  thin 
oxide  layer  in  the  N-S  interface  plays  a  very  important 
role.  Also  experiments  by  Van  Gurp  (1963)  on  sandwiches 
containing  aluminum,  in  which  no  proximity  effect  was 
observed,  may  be  explained  by  the  presence  of  a  thin 
oxide  layer  rather  than  deformed  structure  of  metals 

at  the  interface.  Experiments  by  Miles  and  Smith  (1963) 
on  Al-In  sandwiches  (in  which  Tc  increased  as  the 
sandwiches  were  aged)  may  be  explained  on  the  basis  of 
decoupling  the  metal  films  by  the  introduction  of  a 
gas  layer  in  the  interface. 

(c)  As  pointed  out  first  by  Hilsch  (1962)  ,  the  resulting 
critical  temperature  of  a  proximity  sandwich  is  strongly 
dependent  on  the  mean  free  path  in  the  normal  metal. 

The  mean  free  path  is  also  embedded  in  the  DeGennes- 
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Guyon-Werthamer  theory  for  dirty  superimposed  metals.  In 
McMillan's  theory  the  exact  theoretical  dependence  of 
proximity  effect  on  mean  free  path  is  not  known,  but 
for  constant  ratio  of  the  mean  free  path  to  the  film 
thickness  the  effect  is  roughly  insensitive  to  mean 
free  path.  The  mean  free  path  is  usually  measured 
from  the  residual  low  -  temperature  resistivity,  but 
for  thin  evaporated  films  a  few  assumptions  have  to  be 
made  which  are  difficult  to  verify. 
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CHAPTER  III 

THE  MCMILLAN  TUNNELING  MODEL  OF  PROXIMITY  SANDWICHES 

III.l  Description  of  the  Model 

The  following  assumptions  are  made  about  the  experimental 
system : 

(a)  A  film  of  superconductor  S  of  thickness  dg  is  separated 
by  a  potential  barrier  with  electron  transmission 
probability  a  from  a  film  of  normal  metal  N  (or 
different  superconductor)  with  thickness  d^.  The 
electron  transmission  probability  a  includes  not  only 
the  effect  of  the  oxide  layer  but  also  that  of  any 
impurity  or  imperfection  at  the  interface  of  the  two 
films.  Thus  even  a  "clean"  interface  might  have  a  <1. 

(b)  The  tunneling  Hamiltonian  is  used  to  describe  the 
penetration  of  electrons  though  the  barrier,  which 
restricts  the  transmission  probability  a  to  be  much 
less  than  one. 

(c)  Both  N  and  S  are  thin  compared  with  the  coherence 
length. 

(d)  The  ratio  of  the  mean  free  path  i  to  the  film 
thickness  is  greater  than  or  equal  to  unity. 

Using  this  model,  McMillan  has  derived  the  simultaneous 

T- 

self-energy  equations  for  fixed  BCS  potentials 
ANph  (in  films  S  and  N  respectively)  to  be 


and 
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VE)  = 


As(E)  = 


A 


ir 


N 


[E2 - A2 (E) ] 

ir„ 


1/2 


III .  1 


[E2-A2 (E) ] 1/2 


The  self-consistency  conditions  for  the  BCS  potentials  at 
temperature  T  are 


ph 


ph 


B 


X  /  Re{  - 

N  A  [E2-A2 (E) ] 


an(e)  1  ,  E 

172  }tanh  2kT  dE' 


N 


D 


A  (E) 

X  /  Re{  - - p-r- 

b  C  [E  2  -  A2  (E)  ]  1/2 


III.  2 


}  tanh 


2kT 


dE 


The  correct  integration  limits  are  Ap^  and  [w^,+ ( A^^1)  2  ] 
but  the  contribution  to  the  integrals  is  zero  between  zero 


and 

Aph, 

and 

Aph 

<<  l so  the  limits 

may  be  written  0  and 

WC ' 

thus 

B  = 

N 

5 

0  =  0)^,  and  A  =  C  = 

0. 

The  electronic  density  of  states  can  be  computed  from 


ns,n  =  Re  {  |E|/!E2-A|(N(E)]1/2}. 


III. 3 


In  the  above  equations  F^  and  F^  fulfil  the  conditions 
correctly  stated  by  Kaiser  and  Zuckermann 


Tn  =  TTM2AdsNs(0)  =  h/2xN, 

Ts  =  TTM2AdNNN(0)  =  h/2xs, 

where  M  is  the  tunneling  matrix  element,  A  is  the  area  of 
the  junction,  Nv  ,-.(0)  is  the  density  of  states  at  the  Fermi 
level  in  the  normal  state,  and 


are  relaxation  times 


**  •  —  *. 
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in  the  normal  and  superconducting  material  respectively. 
The  average  time  t  which  an  electron  spends  in  the 

IN  f  o 

N  or  S  film  respectively  may  be  written  as 


TN,S  LN,s/VFN,Sa  ' 


where  v„-T  c  is  the  Fermi  velocity  and  L_T  is  the  average 

r  N  r  o  N  r  b 

electron  path  length  between  collisions  with  the  barrier 
in  the  N  or  S  film  respectively. 

L  =  4d  c  if  the  following  conditions  are  fulfilled: 
Cl)  the  mean  free  path  £  is  infinite  and  (2)  the 

1 NJ  f  O 

probability  that  an  electron  penetrates  the  barrier  to 
enter  the  film  at  an  angle  0  to  the  normal  obeys  the 
cosine  law  (see  Appendix  I) .  For  a  film  which  does  not 
fulfil  these  conditions,  we  may  write 


Lm  c  =  2B 
N,  S 


(  )d 

N's  Vs  N's  ' 


where  B  is  a  function  of  the  ratio  of  the  mean  free  path 
to  the  film  thickness.  Also,  we  can  write 


=  fivFNa 
N  “  4BNdN 


and 


=  ^VFSa 

s  4Bsds 


III .  4 


A  maximum  practical  value  for  these  parameters  is 
T  k  10  meV.  This  is  obtained  with  d  -  15  8  (where  our 

films  became  continuous) ,v  -  108  cm/sec,  BN  =  2 ,  and 

a  -  0.1. 

When  it  is  assumed  that  vpN  =  vps  and  B^  =  Bg ,  then 
rS/rN  =  dN/dS  and  rN,S  “  a/dN,S‘ 


" 
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These  simplified  relations  facilitate  a  physical 
understanding  of  the  results  of  computations  using  the 
McMillan  model. 

III. 2  Computer  Calculations  with  the  Model 
for  Aluminum  -  Tin  Sandwiches 

We  have  carried  out  detailed  computer  calculations 
for  the  aluminum-tin  proximity  sandwich  with  an  analysis 
of  the  tunneling  density  of  states  from  both  sides  of  the 
sandwich.  The  results  are  presented  here  in  graphical  form. 
We  refer  to  aluminum  as  the  normal  metal  (even  if  it  is 
a  superconductor)  and  tin  is  the  superconductor.  The 
calculations  have  been  done  for  the  thickness  range  d  <  dc 
for  the  aluminum  side  of  the  sandwich  and  for  the  thickness 
range  d  >  d  for  the  tin  side  of  the  sandwich. 

IN  o 

The  interesting  feature  of  the  McMillan  model  is  the 
double-peaked  density  of  states  that  appears  in  the  aluminum 
side  of  the  sandwich  when  d^  <  dg  (Tg/r^  <  ^  anc^  a 
small,  and  in  the  tin  side  of  the  sandwich  for  all  ratios 
of  the  thicknesses  as  soon  as  a  is  sufficiently  small.  In 
the  N-side  when  d^  approaches  dg  (Tg/r^  ->  1)  the  peak  P2 
(see  Fig. III. 1)  is  very  small  and  not  observable,  but  when 
a  is  increased,  the  density  of  states  exhibits  a  broad 
single  peak.  The  exact  shape  and  energy  position  of  the 
peaks  depends  on  parameters  Tg,  Ag*1 ,  and  .  The 

situation  in  the  N-film  is  illustrated  in  Fig. III. 1  for 
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FIG.  III.l 


The  McMillan  normalized  density  of  states  in  the 

N-side  of  an  aluminum  -  tin  proximity  sandwich  at  fixed 

values  of  the  potentials  =  l.OmeV,  =  0.1  meV  and 

8  -1 

B  =  B =  2 ,  v„,,  =  =  10  cm  sec  and 

N  S  FN  FS 
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For  higher  a  the  peaks  are  closer  together  and  larger  dg 
is  required  to  resolve  the  double- peaked  structure. 
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two  different  values  of  a,  and  for  the  S-film  in  Fig. III. 2 
for  three  different  values  of  a.  The  curves  in  Figs. III. 1 
and  III. 2  were  calculated  using  eqs.III.l,  fixed  potentials 
AS  ~  ^  me^ '  and  =  0.1  meV ,  and  the  simplifying 

conditions  Bg  =  BN  =  2  and  vpN  =  vps  =  108  cm/sec. 

The  appearance  of  the  double-peaked  structure  may  be 
understood  qualitatively  in  terms  of  two  effects  that  the 
presence  of  each  metal  has  on  the  other.  First,  the  wave 
functions  of  the  excitations  in  the  tin  are  coupled  into 
the  aluminum  creating  a  peak  in  the  aluminum  density  of 
states  near  A^*1,  and  the  wave  functions  of  the  excitations 
in  aluminum  are  coupled  into  tin  creating  a  peak  in  the 
tin  density  of  states  near  APh;  second,  as  a  increases 
the  pairing  energies  in  the  two  metals  move  toward  a 
common  value  determined  largely  by  the  thicker  metal. 

Thus  for  dg  >  d^  the  low  energy  excitation  peak  in  aluminum 
moves  to  higher  energies  as  ct  increases,  and  also  the 
eigenstates  of  the  system  become  less  well  defined,  so  that 
the  peaks  gradually  broaden  and  merge.  Similarly  for 
d^  >  dp  the  high  energy  excitation  peak  in  tin  moves  to 
lower  energies  as  a  increases. 

The  temperature  variation  of  the  peak  positions  and  of 
the  BCS  potentials ,  and  A^,  can  be  calculated  by 
combining  the  self-energy  eq.III.l  with  the  self-consistent 
condition  expressed  by  eq.III.2.  For  aluminum  and  tin  the 
bulk  values  A^  =  0.171,  =  32.21  meV ,  Ag  =  0.246,  and 

g 

ojc  =  16.78  meV  were  substituted.  The  situation  for  dN  <  dg 
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FIG.  III. 2 

The  McMillan  normalized  density  of  states  in  the 
S-side  of  an  aluminum  -  tin  proximity  sandwich  at  fixed 
values  of  the  potentials  A  ^  =  1.0  meV ,  A  ^  =  0.1  meV 

o  IN 

ft  -i 

and  B.7  =  B0  =  2 ,  v_,.T  =  =  10  cm0  sec  and 
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V8) 

rg  (meV) 

rN(meV) 

rS/rN 

a 

a 

75 

7.5 

0.02 

0.2 

0.1 

0 . 002 

b 

75 

75 

0.02 

0.02 

1 

0.002 

c 

75 

750 

0.02 

0.002 

10 

0.002 

d 

75 

75 

0.2 

0.2 

1 

0.02 

e 

75 

750 

0.2 

0.02 

10 

0.02 

f 

75 

7500 

0.2 

0.002 

100 

0.02 

g 

75 

750 

2.0 

0.2 

10 

0.2 

h 

75 

7500 

2.0 

0.02 

100 

0.2 

At  a  given  value  of  a  the  shape  of  peak  structure  is  almost 
insensitive  to  the  actual  thickness  of  the  films.  For 
clarity  the  origin  N  =  O  is  displaced  for  different  curves. 
The  results  of  this  figure  and  Fig. III. 1  are  not  to  be  mixed 
with  the  self consistent  discussion  in  page  20. 
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is  shown  in  Fig. III. 3.  The  calculations  were  done  for  the 
ratio  r  g/rN  =  0.1.  It  can  be  seen  that  for  weak  coupling 
between  the  films  (case  a)  the  higher  energy  peak  P2  ~ 

O 

follows  closely  the  BCS  temperature  variation  of  the  tin 
gap,  while  the  lower  energy  peak  PI  and  at  low 

temperature  follow  the  BCS  temperature  dependence  for  the 
aluminum  gap,  but  deviate  progressively  as  the  temperature 
is  raised  and  finally  approach  zero  at  of  the  overall 
sandwich.  As  the  interaction  between  the  films  is  made 
stronger,  (cases  b  and  c  in  that  order) ,  the  value  of 
P2  -  is  lowered  for  all  temperatures,  while  PI  and 

A^1  are  shifted  to  higher  energies.  The  critical 
temperature  of  the  sandwich  is  very  close  to  the  critical 
temperature  of  pure  tin. 

The  temperature  variation  of  the  peak  positions,  of 
A^h,  and  A^h  for  d  >  d  is  shown  in  Fig. III. 4.  The 
calculations  were  done  for  the  ratio  rg/rN  =  5. 

Qualitatively  the  same  general  behaviour  of  the  peak 
positions,  A^h,  and  A^h  as  in  Fig. III. 3  is  observed,  but 
in  this  case  the  critical  temperature  of  the  sandwich  is  very 
much  different  from  the  critical  temperature  of  pure  tin. 

When  the  superconductor  is  very  little  disturbed  by 
proximity  effect  (that  is,  either  d^  <<  dg  or  a  is  very 
small)  the  temperature  dependence  of  A^h  is  well 
approximated  by  the  BCS  relationship  for  A(T)/A(T=0) , 
so  that  the  critical  temperature  Tc  of  the  sandwich  may  be 
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FIG.  III.  3 


Temperature  variation  of  the  peaks  PI  and  P2  in 

the  density  of  states  in  the  aluminum  side  of  an  aluminum- 

tin  proximity  sandwich  and  the  potentials  AP*1  and  AP*1 

S  N 

when  <  dq . 

Bulk  values  =  32.21  meV ,  A  =  0.171  were  used 

c  N 

for  aluminum;  ^  =  16.78  meV ,  A  =  0.246  for  tin. 

C 

rg/rN  =  0.1  for  all  curves,  =  0.005  meV  for  a, 

F  =  0.05  meV  for  b,  ^  =  o.5  meV  for  c.  AP^1  is  usually 

indistinguishable  from  P2 .  BCS  variation  of  the  tin  gap 

is  indistinguishable  from  A^.  T  =  3.80°K  for  a, 

o  a  c 

T  =  3 . 77 °K  for  b,  T  =  3.69°K  for  c.  These  critical 

k/  v— • 

temperatures  were  obtained  using  eq.III.5. 

-  BCS  variation  of  the  aluminum  gap. 
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FIG.  III.  4 


Temperature  variation  of  the  peaks  PI  and  P2  in 
the  density  of  states  in  the  tin  side  of  an  aluminum  - 
tin  proximity  sandwich  and  the  potentials  and  A?*1 

when  dN  >  dg. 

Bulk  values  =  32.21  meV ,  A  =  171  were  used 

c  N 

for  alumi  num:  oo  =16.78  meV .  A_  =  0.246  for  tin. 


c  '  b 

r^/T  =  5  for  all  curves,  =  0.005  meV  for  a, 

Tg  =  0.05  meV  for  b,  and  V g  =  0.5  meV  for  c.  BCS 
variation  of  the  tin  gap  is  indistinguishable  from 
Tc  =  3 . 80 °K  for  a  ,  Tc  =  3.42°K  for  b  ,  and  Tc  = 
for  c  . 


1. 78°K 


BCS  variation  of  aluminum  gap. 


- 


BCS 


'O  ro 
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found  from  the  BCS  relationship 

2 Agh (T  =  0)  /  kTc  ~  3.5  III. 5 

This  is  the  case  for  the  curves  in  Fig. III. 3  and  for  curve 
a  in  Fig. III. 4.  When  the  conditions  for  the  approximation 
III. 5  are  not  fulfilled,  we  may  estimate  directly  from 
plots  similar  to  Figs. III. 3  and  III. 4,  or  use  the 
approximation  suggested  by  Mr.R.Teshima  (see  Appendix  II)  / 
which  enables  us  to  calculate  T  directly.  (Direct 
calculation  of  from  the  self-consistent  eqs.III.2  cannot  be 
done  because  convergence  of  the  iteration  method  that  was 
used  is  very  slow  around  T  ).  The  results  of  Teshima's 
approximation  (full  line)  and  the  estimation  of  T^,  from 
Figs. III. 3  and  III. 4  or  from  the  approximation  III. 5,  when 
applicable,  are  plotted  vs  T N  in  Fig. III. 5.  One  can  see 
that  Teshima's  approximation  and  the  approximation  III. 5 
agrees  within  2%  (curves  for  r C/T  .  =  0.1,  0.05,  and  0.01). 
Estimation  of  T^  from  plots  similar  to  Figs. III. 3  and  III. 4 
differs  from  Teshima's  approximation  at  most  by  about 
0.2°K  (curves  for  r g/r n  =  0-5,  1,  5,  and  10) ,  but  one  has 
to  keep  in  mind  the  poor  accuracy  of  the  Tc  estimated  from 
Figs. III. 3  and  III.  4.  We  regard  Teshima's  approximation  as 
satisfactory . 

Each  line  of  constant  F g/r N  in  Fig. III. 5  is  a  line  for 
constant  ratio  of  thicknesses  d^/d^.  may  then  be 
increased,  while  keeping  d constant,  by  the  increasing  a  , 
that  is,  producing  a  cleaner  interface.  We  see  that  at 
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FIG.  III. 5 


The  critical  temperature,  Tr  ,  of  an  aluminum-tin 
proximity  sandwich  plotted  vs  T ^  .  The  parameter  is 

rQ/r. 

S  N 


Teshima's  approximation, 

estimate  from  plots  similar 
to  Figs. III.  3  and  III. 4, 

estimate  using  approximation  eq.III.5. 

The  meaning  of  the  symbols  is  as  follows: 


VrN 

Ts/T 

0.01 

A 

1.0 
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fixed  thicknesses,  and  dg,  the  critical  temperature 
approaches  a  constant  limiting  value  as  the  interface  is 
made  cleaner  (higher  r^) .  The  value  of  T  when  the 
interface  is  clean,  therefore  depends  only  on  the  ratio 
of  thicknesses  d^/dg  (-  and  does  not  depend  on  a- 

This  behaviour  of  T^  is  similar  to  the  behaviour  of  the 
Cooper  model  of  the  proximity  effect  with  a  perfectly 
clean  interface  (see  eq.II.2) ,  which  is  also  dependent  on 
the  ratio  of  thicknesses  only. 

When  self-consistent  values  of  and  A^,  and  of  peak 

positions  and  peak  separations  in  the  N  side  of  the 
proximity  sandwich  are  plotted  against  for  different 

ratios  rc/rM,  one  finds  a  very  important  result.  When 
dN  <  dg  the  properties  of  the  system,  shown  in  Fig. III. 6a, 
b,  c,  are  independent  of  (and  therefore  of  Bg  or  dg) 
over  a  wide  range  of  r  .  This  feature  is  used  for 
evaluation  of  our  experimental  data.  In  Fig. III. 6d 
A^h  and  are  shown  as  functions  of  T  for  different 

ratios  of  r q/T m  to  emhasize  that  these  properties  are  much 

O  IN 

more  dependent  on  than  on  F  g . 

On  the  other  hand,  when  dg  <  dN,  the  situation  is 
reversed  and,  as  shown  in  Fig. III. 7a,  A  and  A^  are 
relatively  independent  of  F  N  (  and  therefore  of  or  d^) . 

In  Fig. III. 7b  these  quantities  are  plotted  against  r N  to 
emphasize  that  they  depend  more  strongly  on  than  on 
In  Fig. III. 7c  and  d  the  peak  positions  in  the  S-side  of  the 
proximity  sandwich  are  plotted  vs  F  g  at  1.0°K  and  at  1.4°K. 
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FIG.  III. 6 

The  s elf -consistent  calculations  of  the  McMillan 

model  for  an  aluminum-tin  sandwich  when  d  <  d0 

N  S  • 

(a)  Peak  positions  in  the  N  side  as  a  function  of  r  . 

N 

(b)  Peak  separation  in  the  N  side  determined  from  (a) 

as  a  function  of  F  . 

N 

(c)  Values  of  potentials  a?*1  and  A^n  as  a  function 

°f  rN. 

(d)  the  same  as  (c)  but  as  a  function  of  T  . 

O 

The  properties  of  the  system  are  weakly  dependent 
on  T g  whereas  the  strong  dependence  on  is 

emphasized  in  (d) . 

All  curves  are  plotted  at  1.4°K. 


— •  —  * 
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FIG.  III. 7 

The  self-consistent  calculations  of  the  McMillan 

model  for  an  aluminum  -  tin  sandwich  when  dN  >  dg. 

(a)  Values  of  potentials  and  as  a  function 

of  at  1 . 0°K. 

(b)  The  same  as  (a)  but  as  a  function  of  T  . 

(c)  Peak  positions  in  the  S  -  side  as  a  function  of 
rs  at  1 . 0 °K . 

(d)  Peak  positions  in  the  S  -  side  as  a  function  of 
rg  at  1 . 4 °K . 

^  =  Q  ~  O’  ~  an<^  O  ^g/F^ 
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The  dependences  in  Figs. III. 6  and  III. 7c  are  plotted  for 
a  temperature  of  1.4°K  (above  the  critical  temperature  of 
aluminum)  while  in  Fig. III. 7a,  b,  and  d  they  are  plotted 
for  1°K,  which  is  below  the  critical  temperature  of  aluminum, 
Therefore  and  peak  positions  approach  zero  at  extremely 

weak  coupling  (small  r)  in  Figs. III. 6  and  III. 7c  whereas 
they  approach  the  value  AA1  (T  =  1°K)  =f  0.13  meV  in  Fig. III. 
7a,  b,  and  d. 

Next  it  is  useful  to  consider  the  conditions  under  which 
the  double  peak  in  the  density  of  states  may  be 
experimentally  observed. We  may  deduce  from  Figs. III. 6a  and 
III. 7a  that  two  separated  peaks  exist  for  rN  ~  1  meV  in 
the  N-side  of  the  proximity  sandwich  and  for  rc  5  1  meV  in 
the  S-side  of  the  proximity  sandwich.  For  higher  T's  the 
peaks  merge.  Lower  bounds  for  each  T  may  be  estimated  from 
the  plot  of  peak  magnitudes  in  Fig. III. 8.  If  we  consider 
0.01  change  in  the  normalized  density  of  states  as  the 
minimum  resolvable  change  (because  of  large  variations  of 
the  density  of  states  in  the  gap  region,  smaller  changes 
are  difficult  to  observe) ,  then  we  can  set  the  following 


limits , 


0.005  meV  <  5  1  meV ,  N-side, 


0.001  meV 


r 


1  meV  , 


S-side . 


Using  eqs.III.4  and  vpN  =  vpq  =  10 8 cm/sec  and  = 


FS 


N 


III.  6 


Bs  =  2  ' 


we  may  rewrite  eg. III. 6 


in  the  form 
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FIG.  III. 8 


The  magnitude  of  the  induced  peak  at  1°K  according 
to  the  McMillan  model.  The  left  hand  side  refers  to  peak  P2 
created  in  aluminum  by  the  presence  of  tin  and  the  right 
hand  side  refers  to  peak  PI  created  in  tin  by  the  presence 
of  aluminum.  The  magnitude  of  the  peaks  means  the  distance 
from  the  dip  between  the  two  peaks  to  the  corresponding 
peak  maximum,  as  shown  in  the  inserts. 

Note  that  the  tin  peak  induced  in  aluminum  is 
already  very  small  for  ^g/T  =  1,  while  the  aluminum  peak 
induced  in  tin  is  reasonably  large  for  all  ratios  of  Tg/r^. 

The  scatter  is  caused  by  calculation  of  the  density 
of  states  at  discrete  energy  points  only.  The  meaning  of 
the  symbols  is  the  same  as  in  Figs. III. 6  and  III. 7. 
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6x1065~$  10  3, 
dN 


1  x  10~6  ~  ~  ~  10  3, 

ds 


N-side , 


S-side , 


III.  7 


where  d^  or  dg  are  expressed  in  S.  If  e.g.  d^  -  100  R,  then 

we  should  be  able  to  observe  the  double  peaked  structure  in 

—  4  <  < 

the  N-side  for  a  m  the  range  6x10  ~  a  ~  0.1,  which 

agrees  reasonably  with  our  experimental  observations. 
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CHAPTER  IV 


EXPERIMENTAL  PROCEDURE 


IV . 1 .  Sample  Preparation 

All  tunnel  junctions  were  prepared  on  Fisher  brand 
12  -  550  glass  microscope  slides,  which  were  found  to  have 
very  smooth  surfaces.  The  slides  were  washed  with  detergent 
and  then  slightly  flame  polished.  Contacts  were  soldered  to 
the  glass  slides  with  indium  and  coated  with  gold  prior  to 
the  actual  junction  preparation.  Some  of  the  slides  were 
also  precoated  with  SiO  (Strongin  at  al  (1970) )  in  order  to 
obtain  continuous  metallic  films  at  smaller  thicknesses. 

The  effect  of  SiO  precoating  is  discussed  later.  The  films 
of  metals  were  evaporated  from  a  distance  of  about  50  cm 
through  wedge  shaped  masks  which  were  pressed  against  the 
substrate.  It  was  found,  using  the  Tolansky  interferometer, 
that  wedge  shaped  masks  produced  films  with  more  rectangular 
cross-section  than  "chimney"  shaped  masks.  The  width  of  the 
films  was  ~  0.8  mm.  Aluminum  was  evaporated  from  a  helix  or 
basket  made  from  stranded  tungsten  wire;  tin  was  evaporated 
from  a  tantalum  boat.  The  evaporation  rate  for  tin  was  from 
50  to  100  S/sec , for  aluminum  from  5  to  10  8/sec. 

The  vacuum  system  consisted  of  a  6  inch  oil  diffusion 
pump  with  a  liquid  nitrogen  cold  trap  backed  by  a  mechanical 
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pump.  During  evaporation  liquid  nitrogen  was  also  passed 

through  a  coil  in  the  evaporation  space.  The  vacuum  was 

-7 

about  1  x  10  Torr  before  evaporations  were  started  and 
usually  did  not  exceed  10  Torr  during  the  evaporation. 

IV . 2  Oxidation  (the  Formation  of  a  Tunneling 
Barrier) 

The  standard  procedure  during  the  tunnel  junction 
preparation  was  as  follows.  The  base  layer  of  aluminum 
was  evaporated  onto  the  glass  substrate  and  oxidized  to 
form  a  tunnel  barrier;  then  cross  strips  of  another  metal 
were  evaporated  over  the  oxidized  aluminum  to  form  tunnel 
junctions . 

To  oxidize  the  aluminum  base  layer  we  chose  a  glow 
discharge  oxidation  method.  A  quantitative  discussion  of 
glow  discharge  oxidation  of  evaporated  aluminum  films  for 
a  particular  geometry  of  glow  electrode  is  presented  in 
Appendix  III.  For  different  geometries  of  the  glow  electrode 
the  oxidation  rates  were  different,  but  the  oxidation  laws 
were  the  same.  In  most  cases  the  glow  electrode  was  a  ring  of 
aluminum  wire  2  inches  in  diameter ,  located  about  3  inches 
from  the  substrate.  The  aluminum  strip  was  oxidized  for 
about  350  sec  in  an  atmosphere  of  dry  oxygen  at  a  pressure 
of  0.15  Torr.  The  electrode  was  maintained  at  a  negative 
potential  of  about  350  V  and  the  glow  current  was  one  mA. 

This  oxidation  procedure  produced  a  junction  resistance 
about  lOOq  for  an  A1-A12C>3-A1  junction  if  the  aluminum  of 
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the  proximity  sandwich  was  evaporated  within  a  few  minutes 
after  the  oxidation  and  was  at  least  100  8  thick.  Thinner 
layers  required  less  oxidation  time  to  produce  the  same 
resistance  whereas  a  delay  of  about  2  hours  before 
evaporating  the  aluminum  cover  layer  necessitated  twice  the 
oxidation  time.  If  tin  was  evaporated  onto  the  aluminum  oxide, 
the  reduction  of  the  oxidation  time  was  found  to  agree  with 
Handy's  (1962)  measurements. 

The  actual  oxidation  time  depended  on  such  factors 
as  the  evaporation  rate  of  the  aluminum,  vacuum  during  the 
evaporation,  aluminum  thickness,  the  kind  of  substrate 
(glass  or  SiO)  and  temperature.  The  exact  relationships 
between  these  effects  and  the  oxidation  time  necessary  to 
produce  a  good  junction  were  not  studied  but  an  attempt 
was  made  to  keep  them  the  same  for  every  junction  prepared. 
Then  the  oxidation  by  glow  discharge  was  very  reproducible 
and  the  junction  resistivity  could  be  predicted  with  about 
20%  accuracy. 

The  usual  method  for  estimating  the  oxide 

thickness  during  the  oxidation  is  to  measure  the  oxidation 
time  and  to  use  the  assumption  that  the  resistivity 
increases  linearly  with  time  (Appendix  III.).  The 
measurement  of  the  resistance  change  of  the  aluminum  film 
which  is  being  oxidized  would  provide  an  alternate  method 
for  determinating  the  oxidation  progress.  Assume  that  a 
slab  of  material  has  before  oxidation  a  thickness  d  and 
resistance  R  (see  Fig. IV. la).  During  oxidation  a  thickness 
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FIG.  IV. 1 

The  relative  change  of  the  resistance  of  an  aluminum 
film  during  the  glow  discharge  oxidation  as  a  function  of 
thickness . 

(a)  Configuration  of  the  metal  film  and  oxide  used 
to  estimate  the  resistance  change  during 
oxidation. 

A  R 

(b)  O  experimental  —  vs  d, 

-  plot  of  eq.  IV . 1  for  5  =  3  8  / 

plot  of  eq.  IV . 1  for  6  =  6  S. 
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6  of  aluminum  metal  is  used  to  build  the  oxide  and  the 
resistance  increases  by  AR.  Then,  assuming,  that  the 
electron  scattering  at  the  oxidized  surface  is  unchanged  by 
the  oxidation,  we  may  write 


AR/R  5_  or  AR  =  6_ 
1+A  R/R  d  R  d 


IV.  1 


if  AR/R  <<  1  (which  is  usually  the  case).  A  plot  of  eq.IV.l 

for  6=  3  8  and  <$ =  6  S  is  compared  with  experiment  in  Fig. 

IV. lb.  For  thicknesses  below  200  S  the  oxidized  thickness 

of  aluminum  is  6  -  3  8,  for  higher  thicknesses  6  -  6  S. 

On  the  average  4  to  5  R  of  aluminum  metal  are  used  for 

oxidation.  The  scatter  in  the  experimental  results  may  be 

caused  by  many  effects  -  e.g.  change  of  strain  in  the  film, 

annealing  during  the  oxidation,  or  change  in  the  kind  of 

electron  scattering  at  the  surface  as  it  is  oxidized.  When 

5  S  of  aluminum  is  completely  oxidized,  then  using  the 

oxide  density  p  =3.17  g/cm  (Dignam  et  al  (1966))  and 

assuming  perfect  Al^O^,  one  obtains  ^9  S  for  the  oxide 

thickness.  A  direct  measurement  of  the  mass  gain  of  oxygen 

during  the  oxidation  (see  Appendix  III)  gives  %27  2  for 

the  oxide  thickness  (again  assuming  perfect  Al^O^  and 

o  =3.17  g/cm3).  The  difference  indicates  that  the 
M  ox 

composition  of  oxide  is  not  A^O^.  If  we  assume  that  5  S 
of  aluminum  metal  (0.135  y  g/cm  )  is  combined  with  0.4  yg/cm 
of  oxygen  to  form  an  oxide  (corresponding  to  27  R  of  oxide 
with  density  3.17  g/cm3)/  then  instead  of  having 


. 
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stoichiometric  composition  Al^O^  we  shall  get  composition 
or  <~>ne  ^as  rememt)er  inaccuracies  of  the 

methods;  e.g.,  annealing  during  the  oxidation  should  decrease 
the  resistance  of  the  aluminum  film,  therefore  the  resulting 
resistance  change  is  smaller  than  without  annealing  and  the 
calculated  content  of  oxygen  is  greater. 

IV. 3  Thickness  Measurements 


IV . 3 . a  The  Tolansky  interferometer 

The  basic  tool  for  the  measurement  of  the  thickness 
of  thin  films  is  the  Tolansky  multiple  beam  interferometer, 
described  by  Tolansky  (1948),  Tolansky  (1960)  and  Heavens 
(1951) .  When  the  thin  wedge  of  air  formed  by  two  reflecting 
mirrors  (one  of  them  partially  transparent)  is  placed  in 
the  path  of  a  light  beam,  the  interference  fringes  arising 
from  path  differences  between  the  light  reflected  at  the 
two  surfaces  can  be  observed.  A  narrow  thin  film  evaporated 
onto  the  surface  of  the  opaque  mirror  produces  a  shift  of 
the  position  of  the  fringes  as  shown  in  Fig. IV. 2.  The  height 
of  the  step  A  produced  on  each  fringe  at  the  edge  of  the 
thin  film  is  related  to  its  thickness  d  by  the  following 
expression : 


d 


X  t±_ 

2  Z' 


IV.  2 


where  l  is  the  separation  of  the  fringes  and  X  is  the  wavelength 
of  the  light.  The  actual  experimental  procedure  is  as 
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FIG.  IV. 2 

Illustration  of  the  thickness  measurement  . 

(a)  Cross-section  through  the  interference  element  , 

(b)  The  resulting  fringe  pattern  . 
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follows.  Two  optical  flats  are  cleaned  by  washing  in 
and  then  exposing  them  to  a  glow  discharge . The  film  whose 
thickness  is  to  be  measured  is  evaporated  onto  the  surface 
of  one  optical  flat  and  is  subsequently  coated  by  an  opaque 
layer  of  evaporated  silver.  Another  optical  flat  is  coated 
with  a  layer  of  silver  with  about  94%  reflectivity 
(resulting  in  about  60  or  more  beam  interferences) .  The 
optical  flats  are  placed  together  to  produce  a  thin  wedge 
of  air  as  shown  in  Fig. IV. 2. a. 

Fig. IV. 3  is  a  diagram  of  our  interferometer.  The 
light  from  the  mercury  lamp  is  passed  through  a  green  filter 
and  condensed  onto  an  aperture  of  ^  1  mm  diameter.  A 
parallel  beam  of  light  is  formed  by  a  lens  with  a  focal 
length  of  10  cm.  The  light  then  passes  through  a  half¬ 
transparent  mirror,  interferes  in  the  thin  wedge  of  air 
and  is  reflected  into  the  lens  system  and  camera. 
Manipulation  of  the  thin  wedge  of  air  is  facilitated  by  the 
holder  described  by  Heavens  (1950) .  Photographs  of  the 
fringes  are  evaluated  by  the  method  described  in  Appendix  IV . 

IV . 3 . b  The  Quartz  crystal  thickness  monitor 

The  thickness  of  the  evaporated  films  was  monitored 
during  the  evaporation  by  the  quartz  crystal  monitor.  Two 
5  MHz  crystals  were  mounted  in  the  same  water-cooled 
holder.  One  crystal  was  exposed  directly  to  the  evaporating 
material,  the  other  was  simply  in  thermal  contact  with 
copper  holder  (see  Fig. IV. 4),  resulting  in  at  least  partial 
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FIG.  IV. 3 


Schematic  diagram  of  Tolansky  interferometer  . 
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FIG.  IV. 4 

Quartz  crystal  holder.  The  springs  used  are 
as  weak  as  possible. 
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compensation  for  the  heat  shock  during  evaporation.  The 
heat  shock  produced  a  slow  drift  of  the  frequency  toward 
a  steady  value,  after  the  evaporation  was  completed.  The 
magnitude  of  the  drift  was  proportional  to  the  evaporation 
rate  and  to  the  amount  of  evaporated  material,  factors 
that  are  connected  with  the  amount  of  heat  supplied,  and 
was  zero  when  very  thin  films  were  evaporated.  With  the  two 
crystal  arrangement  the  frequency  drift  was  usually  about 
1%  of  the  total  frequency  change. 

The  circuit  diagram  of  the  quartz  crystal  monitor 
is  shown  in  Appendix  V. 

The  crystal  thickness  monitor  was  placed  about 
32  cm  above  the  evaporation  sources,  while  the  sample  was 
placed  about  50  cm  above  the  sources.  Taking  into  account 
the  geometry  following  Holland  (1963),  the  sensitivity  of 
the  crystal  monitor,  as  calibrated  by  the  Tolansky 
interferometer  was 

^sample  =  0>366  £/Hz  for  A1 / 

-?.ai^P.1.^  =  0.115  2/Hz  for  Sn, 


where  d  ,  is  the  thickness  of  the  film  at  the  sample 
sample 

position  and  Af  is  the  frequency  change  of  the  quartz 

crystal  monitor.  When  the  ratio  of  frequency  change  per  2 

for  Al  and  Sn  was  compared  with  the  ratio  of  the  densities 

n  ,  and  of  Al  and  Sn  it  is  found  that 

M  Al  H  Sn 


. 
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Af 


Al 


Af 


=  0.315, 


Al 


=  0.370. 


Snyi^  ^Sn 


The  difference  between  these  values  may  be  caused  by  a 
different  porosity  of  the  evaporated  films  or  by  a  different 
sticking  probability  for  impinging  metal  atoms  on  the 
optical  flat  or  quartz  crystal  surface. 

People  usually  use  silver-  or  gold-coated  crystals 
(e.g.  Riegert  (1968)).  However,  we  have  only  used  crystals 
coated  with  silver  because  their  frequency  stability  was 
superior  to  those  coated  with  gold. 

The  crystals  usually  stopped  oscillating  after  they 
were  loaded  with  mass  equivalent  to  about  105Hz. 


IV . 4  Description  of  the  Apparatus  and  the  Production 
of  Low  Temperatures 

The  tunnel  junctions  were  prepared  in  the  metal 
evaporator  -  cryostat  apparatus  and  cooled  down  without 
breaking  the  vacuum.  This  was  necessary  because  the  transfer 
of  a  junction  into  a  cryostat  through  the  air  always 
destroyed  them  when  very  thin  metal  films  were  involved.  Also 
when  proximity  sandwiches  were  involved,  reproducible 
results  could  not  be  obtained  if  they  were  exposed  to  air. 
Probably  the  interface  between  the  metals  in  proximity  was 
influenced  by  the  air. 

Fig. IV. 5  is  a  schematic  diagram  of  the  apparatus. 

The  crvogenic  part  was  manufactured  by  Sulfrian  Cryogenics, 


. 
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FIG.  IV. 5 


Schematic  diagram  of  the  evaporator-cryostat. 

A  -  port  for  liquid  He  input  into  the  He  pot 
B  -  outlet  for  the  measurement  of  He  vapor  pressure  in 
the  He  pot 

C  -  connection  to  the  pump 

D  -  port  for  liquid  He  input  into  the  main  He  vessel 
E  -  liquid  helium  level  detector 
F  -  radiation  baffle 

G  -  carbon  resistor  for  liquid  helium  level  sensing 
H  -  guide  tube  for  steel  piano  wire 
I  -  pulleys  for  mask  operation  and  pressing  system 
J  -  ports  for  wires 
K  -  germanium  thermometer 
L  -  sample 

M  -  spring-loaded  balls  for  easy  location  of  mask  position 
N  -  movable  mask 
0  -  pressing  mechanismus 

P  -  movable  shutter  on  LN  temperature  shield 
Q  -  window  for  observation  of  glow  discharge  electrode 
R  -  glow  discharge  electrode 
S  -  liquid  nitrogen  coiled  tube 

T  -  optical  flat  for  calibration  of  quartz  crystal  thickness 
monitor 

U  -  quartz  crystal  thickness  monitor 


■ 
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V  -  window  with  mirror  for  observation  of  evaporation 

sources 

W  -  high  vacuum  pump 

Y  -  shutter 

Z  -  multi-f ilament  carousel 


Z^o  3  ^ 
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Inc.  The  liquid  level  was  detected  by  a  capacitance  level 
detector  (see  Appendix  VI) ,  which  was  immersed  in  either 
liquid  nitrogen  or  liquid  helium.  The  200  cm3  He^  pot  was 
attached  to  the  cryogenic  part  by  a  single  stainless  steel 
tube.  The  helium  and  nitrogen  temperature  shields  were 
connected  to  the  corresponding  baths  by  a  greased  metal 
flange.  Shutters  that  could  be  operated  from  outside  the 
apparatus  normally  covered  the  apertures  through  which  the 
evaporating  metal  reached  the  glass  substrate.  The  substrate 
was  clamped  to  the  bottom  of  the  helium  pot  and  Apiezon  N 
grease  served  as  a  thermal  contact.  The  mask  and  the 
mechanisms  which  could  press  the  mask  against  the  substrate 
were  operated  from  the  outside  by  steel  piano  wire  (0.010 
inch  diameter) .  The  glow  discharge  electrode  was  placed 
outside  the  nitrogen  temperature  shield.  Boats  or  filaments 
were  mounted  in  the  circular  carrousel,  which  permitted 
succesive  evaporation  from  six  independent  sources.  One  of 
the  sources  was  equipped  with  a  feeding  mechanism  so  that 
additional  metal  could  be  added  to  the  source  during 
evaporation  of  thick  films.  A  window  with  a  mirror  was 
provided  for  observation  of  the  evaporation  sources.  A 
typical  run  is  described  in  Appendix  VII. 

Temperatures  in  the  range  between  1.3°K  and  4.2°K 
were  measured  by  a  germanium  resistor  placed  on  the  bottom 
of  the  pumping  pot.  The  resistor  was  calibrated  against 
vapor  pressure  measurements  and  the  experimental  points 


were  fitted  by  a  computer  to  the  formula  (shown  e.g.  by 
Blakemore  et  al  (1970)) 
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R 


A  T-B  eC/T 


f 


where  A,  B ,  and  C  are  constants.  A  table  of  R  vs  T  was 
compiled  in  10  m°K  steps,  which  is  the  limit  of  resolution 
of  our  data  collected  over  about  four  months  of  experiments. 
In  the  range  between  1.05°K  (lowest  attainable  temperature 
in  our  system)  and  1.3°K  the  vapor  pressure  was  compared 
with  the  1958  He4  Scale  of  temperatures  compiled  by  Dijk 
et  al  (1960).  The  agreement  of  the  sample  temperature  with 
the  temperature  reading  was  checked  at  3.8°K  (the  critical 
temperature  of  tin)  and  by  comparison  of  tunneling 
measurements  with  the  BCS  variation  of  tunneling 
conductance . 

The  cryostat  was  wired  with  No. 40  B&S  gauge  copper 
wires  which  were  thermally  anchored  at  liquid  nitrogen 
temperature,  liquid  helium  temperature  and  on  the  He^  pot 
before  they  were  connected  to  the  sample. 

The  conductivity  measurements  were  done  on  a 
tunneling  conductance  bridge  similar  to  the  one  described 
by  Rogers  (1970). 
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CHAPTER  V 


DISCUSSION  OF  MATERIALS  USED 


V . 1  The  Selection  of  Materials 

We  have  selected  aluminum  for  the  base  layer  because 
it  is  easily  oxidized  to  form  a  tunneling  barrier  and  also 
its  critical  temperature  (1.2°K  for  bulk  material  and  higher 
for  thin  films)  is  convenient;  experiments  can  be  done 
readily  with  the  base  layer  either  normal  or  superconducting . 

Aluminum  and  tin  have  been  selected  for  the  proximity 
sandwich.  Again  aluminum  was  selected  because  of  its 
critical  temperature  and  ease  of  oxidation  which  is 
important  in  the  control  of  the  transmission  probability, 
a  .  Tin  was  selected  because  it  becomes  superconducting 
below  3.8°K,  and  within  our  range  of  measurement  we  could 
then  have  neither,  only  one,  or  both  metals  of  the  proximity 
sandwich  in  its  superconducting  state. 

V. 2  The  Problem  of  Interdiffusion 

As  suggested  first  by  Rose-Innes  and  Serin  (1961) 
interdiffusion  may  often  occur  in  proximity  sandwiches. 
However  according  to  Hansen  (1958)  ,  and  Elliott  (1965)  the 
solid  solubilities  for  the  aluminum  -  tin  system  are  very 
restricted  and  no  intermetallic  compounds  are  known.  Also 
in  the  experiments  reported  here  a  thin  oxide  layer  was 


’ 
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formed  between  the  aluminum  and  tin  which  would  further 
reduce  diffusion  effect.  Perhaps  the  strongest  evidence  that 
our  procedures  produced  proximity  structures  with  distinct 
aluminum  and  tin  films  lies  in  the  tunneling  results,  which 
exhibit  several  characteristics  that  would  not  arise  in 
badly  interdif fused  or  alloyed  structure.  An  alloy,  as 
pointed  out  by  Guyon  et  al  (1966)  ,  should  behave  like  a  BCS 
superconductor . 


V . 3  Mean  Free  Path  Considerations 

The  important  condition  for  the  McMillan  tunneling 
model  of  the  proximity  effect  to  be  valid  is  that  the 
electron  mean  free  path  be  comparable  with  the  thickness 


-L  uiit;  j — l  xiTi  3  m  question .  The  me  a  sur  eme  n  t  of  the 

mean  free  path  from  the  residual  resistivity  of  a  very  thin 
film  is  usually  complicated  by  the  presence  of  surface 
scattering.  If  we  consider  pure  bulk  material  with  its 
bgilk  mean  free  path  a  ^  and  bulk  resistivity  p^,  then  the 
following  expression,  shown  e.g.  by  Kahan  et  al  (1960),  is 
valid : 


.  J - -  “1  “1  4  ^  -C  A  1 


p0£0  =  n2/'3'e7  (I?)  “constant,  V.X 

where  n  is  the  number  of  conduction  electrons  per  cm3,  e  is 
the  electron  charge  and  h  is  the  Planck  constant.  If  we 
assume  that  a  thin  layer  is  cut  from  the  bulk  piece  of 
material,  without  otherwise  changing  the  structure,  then 
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the  specific  resistivity  of  the  thin  layer  will  be  larger 
than  the  bulk  value  due  to  the  additional  scattering  at  the 
surface.  The  resistivity  increase  as  a  function  of  thickness 
was  first  calculated  by  Fuchs  (1938)  and  recalculated  by 
Sondheimer  (1952).  The  experimental  results  of,  e.g. 

Andrew  (1949)  on  very  pure  rolled  thin  tin  foils  and  on 
mercury  wires  agreed  quite  well  with  the  predictions  of  the 
theory.  A  pure  material  has  very  long  mean  free  path  and 
therefore  relatively  thick  foil  -  which  has  the  structure 
of  bulk  material  -  is  already  affected  by  the  surface 

scattering.  However  the  experiments  of  Hilsch  (1962)  and  of 

S 

Basewitz  and  Minnigerode  (1964)  on  thin  evaporated  films 
cannot  be  fitted  to  the  Sondheimer  calculations.  The 
resistivity  increases  more  rapidly  with  decreasing  thickness 
than  predicted  by  the  calculation.  We  believe  that  in  thin 
films  the  bulk  structure  is  not  preserved  and  as  the  film 
is  made  thinner  more  defects  are  introduced.  Assume  for  a 
moment,  that  we  have  bulk  material  with  the  same  structure 
as  a  thin  film.  Then  its  mean  free  path  is  and  bulk 
resistivity  is  p  .  Further  we  assume,  that  eq.V.l  is  still 

valid  for  the  material  with  the  thin  film  structure  and 
that  (as  suggested  by  Mayadas  and  Shatzkes  (1970)) 

P0l0  =  PbV  V-2 

The  specific  resistivity  of  our  thin  film  is  not  pg  but  P 
(because  of  additional  surface  scattering)  and  we  assume 
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that  p  is  connected  with  p  by  Sondheimer ' s  calculation 

B 

with  the  mean  free  path  Z  . 

B 

Therefore  to  find  the  mean  free  path  Z  of  a  thin 
film  from  its  thickness  d  and  residual  resistivity  p  we 
proceed  as  follows.  The  constant  in  equation  V.l  is  found 
for  a  given  material  (see  e.g.  Bassewitz  and  Minnigerode 
(1964)  who  list  Pq&q  for  various  metals)  and  the  Sondheimer 
dependences  of  p  vs  d  are  calculated  for  different  pairs 
of  pB  and  (which  are  in  accord  with  eq.V.2).  The  family 
of  calculated  p  vs  d  curves  with  as  a  parameter  is 
plotted  together  with  experimental  points  of  p  vs  d  (see 
Fig. V.l)  and  the  corresponding  mean  free  path  is  determined 
as  a  function  of  thickness  d.  For  faster  data  processing 
we  have  written  a  computer  program,  which  for  given  values 
of  p,  d,  and  the  constant  in  eq.V.l,  finds  the  corresponding 
mean  free  path.  All  calculations  have  been  done  for  diffuse 
scattering  at  the  film  surfaces. 

In  Fig.V.2  the  mean  free  paths  for  thin  aluminum 
films  evaporated  at  room  temperature  onto  glass  substrates 
and  glass  substrates  coated  bv  509  R  of  SiO  are  presented. 

1  2  Q 

The  value  pQ£0  =  4.9  x  10  ft  cm 2  was  used  and  the  5  A 
correction  was  applied  to  the  thickness  due  to  the  oxidation 
of  the  surface  (see  Chapter  IV. 2).  It  is  very  interesting  to 
note  that  aluminum  films  evaporated  onto  an  SiO  surface 
possess  a  shorter  mean  free  path  than  the  films  evaporated 
onto  a  glass  surface.  The  glass  substrate  is  usually 
precoated  with  SiO  in  order  to  "improve"  the  quality  of 
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FIG.  V.l 


Residual  resistivity  as  a  function  of  thickness 
for  various  values  of  mean  free  path. 

n  -12  o 

-  calculation  for  A1 ,  P_&  =  4.9  x  10  ficm  , 

o  Jd 

Q  experimental  data,  aluminum  films  evaporated  onto 
glass  substrate  , 

A  experimental  data,  aluminum  films  evaporated  onto 
glass  substrate  coated  by  v  500  £  of  SiO. 
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FIG.  V . 2 


The  mean  free  path  of  evaporated  aluminum  films  as 

a  function  of  thickness  (  evaporation  done  at  room 

12  2 

temperature  ).  The  value  pJD  =  4.9  x  10  ZQ,c rn  and  5  & 

Jd  n> 

thickness  correction  (  due  to  oxidation  )  was  used. 

O  glass  substrate , 

A  glass  substrate  precoated  with  500  R  of  SiO. 

-12  2 

When  the  value  =  8.2  x  10  ftcm  is  used  the  trend 

13  u 

of  vs  d  is  shown  by 

13 

-  for  glass  substrate,  and  by 

for  glass  substrate  coated  with  500  R  of  SiO. 

Solid  lines  represent  visual  fit  to  the  experimental  points. 
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the  films  arid  to  make  them  continuous  at  smaller  thicknesses 

(Strongin  et  al  (1970)).  We  have  actually  observed  that 

films  evaporated  onto  the  SiO  surface  become  continuous 

sooner.  The  effect  is  probably  connected  with  the  formation 

of  droplets  of  evaporated  material  the  size  of  which  depends 

on  the  substrate.  Small  droplets  with  short  mean  free  paths 

are  probably  interconnected  sooner  than  bigger  droplets 

with  longer  mean  free  path.  Slight  oxidation  probably 

maintains  the  partial  separation  of  drops  and  short  mean 

free  path  even  when  the  film  is  thick. 

The  mean  free  path  for  tin  films  evaporated  onto  a 

glass  substrate  is  listed  in  Table  I.  The  value  Pq^q  = 

-12  ? 

10.5  x  10  ficm  was  used.  Note,  that  the  ratio  £  /d  is 

B 

generally  greater  for  tin  than  for  aluminum. 

The  scatter  in  the  data  is  probably  connected  with 

the  fact  that  all  evaporations  were  done  at  room  temperature 

which  results  in  uncontrolled  annealing .Also  the  aluminum  films 

were  oxidized  in  a  glow  discharge  (in  order  to  produce  the 

tunnel  junction)  which  might  change  the  surface  scattering 

and  therefore  the  residual  resistivity. 

The  mean  free  path  fulfils  the  McMillan  assumption 

l  ~  d  even  if  it  is  calculated  by  using  the  smallest  values 
B 

of  p  q£  q  found  in  the  literature.  If  e.g.  the  mean  free  oath  of 

-12  2 

aluminum  is  estimated  using  the  value  =  8.2  x  10  P  cm 

obtained  by  Holwech  and  Jeppesen  (1967)  it  would  be  at 

least  twice  as  big  as  the  mean  free  path  which  we  obtained 

-12  2 

using  the  Chambers'  (1950)  value  pQ£0  =  4.9  x  10  ftcm  . 
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TABLE  I 


The  mean  free  path  of  evaporated  tin  films  as 


function 


of  thickness. 


a  (8) 

V8) 

Vd 

264 

2156 

00 

• 

to 

345 

2382 

6.9 

376 

3345 

8.9 

960 

1570 

1.7 
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CHAPTER  VI 


TUNNELING  INTO  A  PROXIMITY  SANDWICH 


VI . 1  The  Tunneling  Integral 

The  normalized  differential  conductance  a (V)  of  a 
tunnel  junction  formed  between  two  metals  A  and  B  can  be 
expressed  in  terms  of  their  normalized  densities  of 
electronic  states  N..  and  N_.  (assuming  constant  tunneling 

A  n> 

matrix  elements)  by  the  expression 


a  (V) 


(di/dv)s 

(dI/dV)N 


+«>  3N  (E-V) 

/  Na(E)  [ - w -  cj>  (E  , V ) 


+  NR(E-V)^|Ai]dE, 


VI.  1 


where  E  is  energy,  V  is  the  applied  voltage  in  energy  units, 
and  <ME,V)  =  [1+exp  (E-V) /kT]  _1  -  [1+expE/kT]  _1 . 

In  order  to  carry  out  numerical  calculations  the 
functions  NA ,  N„  ,  and  the  derivative  3N  /3V  must  be  finite. 
We  have  used  the  approximation 


where 

with 

of  the 


nb  =  l 


if  B  is  a  normal  metal, 


=  Re [ | E| /  (E2 -A2 ) !/2]  if  B  is  a  superconductor  VI. 2 
A  =  Ab(1  +  i«) 

Ad  the  BCS  energy  gap  in  metal  B. 

N  is  the  McMillan  density  of  states 
A 

proximity  sandwich. 

The  constant  6  has  been  evaluated  by  fitting  eq.VI.l 
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to  experimental  tunneling  data.  6  =  0.013  was  obtained 

by  putting  and  fitting  to  data  for  an 

Al-I-Al  junction  at  T  =  0.4°K,  then  using  this  value  eq.VI.l 
was  again  fitted  to  results  for  an  Al-I-Sn  junction  at 
T  =  1. 0°K,yi elding  6  =  0.038.  A  comparison  of  the 

calculated  and  experimental  curves  for  this  junction  is 
shown  in  Fig.VI.l.a.  The  approximate  density  of  states  for 
a  superconductor  was  then  used  in  eq.VI.l  to  calculate  the 
tunnel  conductance  of  a  normal  metal-superconductor  junction 
Comparison  with  Bermon's  (1964)  calculations  for  a  similar 
junction,  plotted  in  Fig.VI.l.b,  shows  that  the  conductance 
values  agree  within  1%.  We  conclude  that  the  approximate 
density  of  states  of  eq.VI.2  is  adequate  to  use  in  the 
calculation  of  the  tunnel  conductance  of  a  junction  formed 
between  a  superconductor  (or  a  normal  metal)  and  a  proximity 
sandwich  characterized  by  the  McMillan  density  of  states. 


VI . 2  Tunneling  into  the  N-side  of  a  Proximity 
Sandwich 


VI. 2. a  Tunneling  from  a  normal  metal 
A  general  feature  of  tunneling  from  a  normal  metal 
into  the  N-side  of  a  proximity  sandwich  may  be  observed  by 
considering  the  experimentally  measurable  quantity 
[1  -  a  (0)].  For  tunneling  into  a  BCS  superconductor 
[1  -  a  (0)]  a  1  -  T/Tccc  A2  whereas  it  increases  less 


- 

# 


59 


FIG.  VI.  1 


(a)  The  normalized  tunneling  conductance  a (v) 
in  the  junction  Al-I-Sn, 

-  experimental  curve  r 

O  calculation  using  the  smeared  BCS  density  of  states 

(see  eq.VI.2)  with  A^  =  0.22  meV ,  6  ^  =  0.013, 

A  =  0.60  meV ,  and  6  =  0.036. 

bn  Sn 

(b)  The  deviation  of  Bermon's  conductance  a  from 

n> 

the  normalized  tunneling  conductance  a (V)  between  a  normal 
metal  and  the  smeared  superconducting  density  of  states 
(  6=  0.013,  A/kT  =  2.0) . 


" 


- 
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rapidly  than  this  with  decreasing  temperature  just  below 
for  tunneling  into  the  N-side  of  a  proximity  structure 
and  then  grows  faster  as  the  temperature  is  lowered.  Guyon 
et  al  (1966)  proposed  the  empirical  relationship 
[1  -  cr  ( 0 )  ]  a  (1  -  T/T^,) n  which  fits  their  data  for  proximity 
sandwiches  made  with  zinc  and  indium-bismuth  alloys.  They 
and  Claeson  et  al  (1967)  have  obtained  n-values  between  2.3 
and  3.0,  but  their  N-films  were  at  least  1000  S  thick  and 
it  is  probable  that  the  electron  mean  free  path  was  smaller, 
thus  violating  one  of  the  conditions  Of  the  McMillan  model. 

In  Fig. VI. 2  the  behaviour  of  our  aluminum-tin  samples, 
in  which  the  slopes  range  from  about  1  to  6 ,  is  shown. 

Slopes  of  higher  value  are  very  difficult  to  measure  because 
the  temperature  where  the  energy  gap  develops,  which  is 
used  as  the  experimental  ,  cannot  be  measured  accurately 
when  a(0)  is  varying  very  slowly  with  temperature. 

These  results  are  in  good  accord  with  the  McMillan 

model.  For  weak  coupling  (small  a)  the  magnitude  of  is 

very  small,  as  may  be  seen  in  Fig. III. 3,  and  increases 

slowly  with  decreasing  temperature  producing  the  same 

behavior  in  [1  -  a ( 0 ) ]  until  temperatures  near  the  critical 

temperature  of  aluminum  are  reached;  thereafter  these 

quantities  both  rise  rather  sharply.  For  the  strong  coupling 

cases  (case  c.  Fig. III. 3) A  ^  approaches  the  BCS  behavior 

and  so  does  the  experimental  slope  of  [1  -  a(0)].  We 

calculated  a(0)  as  a  function  of  temperature  for  different 

combinations  of  T,,  and  Fc ,  using  the  relationships  of  the 

N  o 
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FIG.  VI. 2 


Experimental  dependence  of  [1  -  a(V=0)]  on 

(1  -  T/T  ) .  a (V=0)  is  the  normalized  tunneling  conductance 
c 

at  zero  bias  between  a  normal  metal  and  the  N-side  of  an 
aluminum-tin  proximity  structure.  The  appropriate 
experimental  values  of  rN  determined  from  Fig.  VI. 3  are  : 
m  rN  =  1.0  meV,  A  rN  =  0.55  meV ,  □  =  0.05  meV ,  and 

0  rN  =  0.02  meV . 


• 

( 0  =  A ) 
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McMillan  model.  Tc  was  estimated  using  eq.III.5  and 

[1  -  cr(0)]  vs  (1  -  T/T^)  was  plotted  logarithmically  as 

shown  in  Fig. VI. 3.  The  curvature  of  these  graphs  is  small 

and  straight  lines  were  readily  fitted  and  their  slopes, n, 

determined.  These  n-values  were  plotted  as  a  function  of 

^°r  ratios  of  r^/T  Just  as  for  the  properties 

shown  in  Fig. III. 6,  the  value  of  the  slope  is  independent  of 

Ts  over  a  wide  range  of  T  .  This  important  result 

facilitates  the  determination  of  I\T  (or  a/B„)  for 

N  N 

experimental  samples  without  the  necessity  of  knowing  r 

u 

(or  B  ) .  The  calculations  with  McMillan's  model  for  r 
b  n 

between  5.0  and  0.001  meV  show  that  n  may  vary  from  1 

(strong  coupling)  to  25  (weak  coupling  ) . 

In  Fig .VI. 4  o(V)  has  been  plotted  as  a  function  of 

applied  voltage  for  sample  D81  (with  d^  =  115  £, 

dg  =  3150  &) ,  together  with  the  curve  for  tunneling  into 

a  BCS  superconductor  and  the  curve  for  the  McMillan  model, 

all  at  T/Tc  =  0.73.  It  is  immediately  clear  that  the 

experimental  results  are  vastly  different  from  the  BCS 

results  and  fit  the  McMillan  model  very  well,  particularly 

when  the  effects  of  a  slight  error  in  the  determination  of 

the  tunnel  conductance  when  the  proximity  sandwich  is  just 

above  Tc  are  considered.  In  order  to  obtain  the  curve  for 

the  McMillan  model,  the  values  =  0.047  meV  and  T  =  2.63°K 

were  used,  and  assuming  B^  =  Bg,  v^N  =  v^  we  obtained 

r  / r  =  0.05.  The  results  are  relatively  insensitive  to  the 
S  N 

value  of  Fc /  hT  t  T  is  the  same  as  that  for  the  experimental 

o  IN 


■  o- 1 1  fef  w  .i-#i  ■  s»i  uiMtfm 
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FIG.  VI. 3 


Calculated  values  of  the  exponent  n  in  the  empirical 

expression  [1  -  g(V=0)]  «  (1  -  T/T  )n  for  tunneling  from 

c 

a  normal  metal  into  the  McMillan  model  for  aluminum-tin 

using  different  values  of  rc/l\T-  In  the  insert  the 

calculated  dependence  of  [1  -  g(V=0)]  on  (1  -  T/T  )  is 

o 

shown. 

1  1  on  m  r-,  1  T  /  r  —  n  i  fa)  r.  =2.0  meV, 

— —  -N/  s  v  '  JN 

(b)  r  =  0.5  meV,  (c)  =  0 . 2  meV ,  and 

(d)  rN  =  0.05  meV. 

-  BCS  dependence. 
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FIG.  VI.  4 


The  normalized  tunneling  conductance  a (V)  between  a 
normal  metal  and  an  aluminum-tin  sandwich. 

0  experiment  D81,  Tc  =  3.6°K,  d^  =  115  8, 

d_  =  3150  8. 

Sn 

-  calculated  from  the  McMillan  model  with 

r  =  0.047  meV,  V  =  0.0024  meV ,  A|h  =  0.474  meV, 

No  b 

and  A^h  =  0.061  meV .  T  =  2.626°K  and  T/Tc  =0.73 
for  both. 

At  zero  bias  only,  the  theoretical  calculations  are 
included  for  r^/r^  =  0.05  and  =  0.04  meV  (  A.  )  , 
rrT  =0.06  meV  (  □  )  . 

For  comparison,  part  of  the  Bermon  conductance  o 
for  the  same  T/Tc  is  shown  by  dashed  line.  crB  (V  =0)  =  0.514. 


. 


a  (V) 
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sample,  and  =  0.047  meV  fits  the  slope  n  =  3.48  (see 
Fig. VI. 3)  which  was  determined  from  a  plot  like  Fig. VI. 2 
for  sample  D81. 

VI. 2. b  Tunneling  from  a  superconductor 
Although  tunneling  into  a  proximity  structure  from  a 
normal  metal  above  1°K  is  clearly  distinguishable  from  that 
into  an  ordinary  superconductor  only  after  some  analysis, 
the  normalized  conductance  for  tunneling  from  a 
superconductor  into  a  weakly-coupled  proximity  sandwich 
exhibits  the  very  characteristic  multiple-peaked  structure 
shown  in  Fig. VI. 5.  b.  This  curve  was  constructed  for 
tunneling  at  1.4°K  between  the  density  of  states  for 
aluminum  and  the  McMillan  density  of  states  shown  in  Fig. 

VI.  5.  a.  The  expected  energies  of  structure  are  indicated  by 
arrows  at  PI  ±  A  and  P2  ±  A.  The  structure  at  P2  -  A  is 
usually  very  small  because  the  concentration  of  thermally 
excited  electrons  which  contribute  most  of  this  peak  and  the 
change  in  the  density  of  states  in  P2  are  usually  much 
smaller  than  in  PI. 

In  Fig. VI. 6  an  experimental  curve  (for  dN  =  100  X, 
d  =  2260  8)  is  compared  with  a  calculated  curve  using 

D 

T  =0.1  meV  and  Ik  =  0.005  meV.  There  are  two  possible 
N  S 

ways  of  fitting  the  calculated  and  experimental  curves. 

From  experiment  the  peak  separation  can  be  measured  and  the 
resulting  obtained  from  Fig. III. 6b,  or  the  position  of 
structure  PI  can  be  measured  (taking  account  of  the 


,  ,  - 
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FIG.  VI. 5 


Illustration  of  the  occurrence  of  multiple  peak 
structure  in  the  tunneling  conductance  between  an  ordinary 
superconductor  and  the  N-side  of  an  aluminum-tin 
proximity  sandwich. 

(a)  The  normalized  densities  of  states  in  the  super¬ 
conducting  aluminum,  (characterized  by  A^  =  0.20  meV , 

6-t  —  u .  0 1  j  i  i  ana  proxiruiLy  sandwich,  N  (characterized 

Al  ^  2  prox 


by  Fg  =  0.025  meV ,  =  0.5  meV , 


Ph  _ 


=  0.56  meV ,  and 


ph 


ANr‘  =  0.31  meV) .  When  voltage  is  applied  the  densities  of 
states  shift  on  the  energy  scale. 

( b)  The  resulting  calculated  tunneling  conductance  at 
1.4°K  (continuous  line).  Expected  structures  at  PI  +  A  and 


P2  +  A  are  indicated  by  arrows. 


For  comparison  the  tunneling  conductance  between 
a  normal  metal  and  the  same  proximity  sandwich  at  the 
same  temperature  is  also  plotted  (dashed  line) . 


-0.5  0.0  0.5 


eV  (meV) 
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FIG.  VI. 6 


Comparison  of  experimental  and  theoretical 
normalized  conductance  for  tunneling  from  superconducting 
aluminum  into  the  N-side  of  aluminum-tin  proximity 
sandwich. 

- —  experiment  E82  ,  T  =  1.4°K,  d^  =  100  X,  and 


d0  =  2260  8. 
Sn 


,ph 


calculated,  F^  =  0.005  meV ,  F^  =  0.1  meV,  -  0.56 

meV,and  A^*1  =  0.22  meV .  Aluminum  electrode  from 
which  the  tunneling  was  done  is  characterised  by 
A^  =  0.20  meV  and  6  =  0.013. 


(V)  N  -SIDE 


V  ( mV ) 
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influence  of  the  gap  of  the  pure  aluminum  electrode)  and 
rN  can  be  determined  from  Fig. III. 6. a.  Structure  PI  is 
selected  because  its  sensitivity  to  is  higher  than  that 
of  P2.  The  first  method  was  used  in  constructing  Fig. VI. 6 
and  values  of  A ^  and  A ^  were  estimated  from  Fig. III. 6. c. 

As  expected  the  fit  of  the  position  of  the  structure  is 
insensitive  to  r  Q  so  that  the  assumptions  B  =  B  and 

o  No 

v  =  v  are  adequate.  We  then  obtain  T  _/T  -  d.T/d„  -  0.0  5 

£  JN  r  o  SNNS 

and  rg  =  0.00  5  meV. 

In  Fig. VI. 7  another  experimental  curve  (for  d^  = 

266  X,  dg  =  2370  S)  is  compared  with  a  calculated  curve 
using  =  1.0  meV  and  Fg  =  0.05  meV .  The  fitting  procedure 
is  the  same  as  for  the  curve  in  Fig. VI. 6.  We  may  see  the 
difference  in  the  shape  of  the  curves  for  different 
transmission  probabilities  a.  For  the  curves  in  Fig. VI. 6 
the  a  is  approximately  10  x  bigger  than  for  the  curves  in 
Fig. VI. 7. 

The  structures  on  the  calculated  curves  are  sharper 

than  the  experimental  structures.  However  our  assumed  density 

of  states  for  aluminum  has  been  shown  to  produce  less  sharply 

peaked  structure  than  experiment  (Fig. VI. 1. a)  so  it  cannot 

be  responsible  for  this  observed  difference.  A  better  fit 

to  the  shape  of  the  structure  in  Fig. VI. 6  would  be  obtained 

using  rs  =  0.05  meV.  This  would  require  BN/BS  -  10  in  order 

to  keep  r  =  0 . 1 .  We  do  not  know  the  exact  theoretical 
c  N 

dependence  of  the  B  values  on  the  ratio  of  the  mean  free 
path  to  the  thickness  ( l/d) .  However  we  may  assume  that  if 
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FIG.  VI. 7 


Comparison  of  experimental  and  theoretical 
normalized  conductance  for  tunneling  from  superconducting 
aluminum  into  the  N-side  of  an  aluminum  tin  proximity 
sandwich . 

- -  experiment  C93  ,  t  =  1.395°K,  d^  =  266  8,  and 

d0  =  2370  R. 

Sn 

-  calculated,  Fg  =  0.05  meV,  =  1.0  meV , 

Agh  =  0.56  meV,  and  =  0.34  meV.  The  aluminum 

electrode  from  which  the  tunneling  was  done  is 
characterized  by  A^  =  0.245  meV  and  6  =  0.013. 


cr ( V )  N-SIDE 


V  ( mV ) 


0 


1.0 
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Vd  0  then  B  goes  monotonically  to  infinity  (B  is 
connected  with  average  path  of  the  electron  in  a  given  film; 
see  App. I) .  Therefore,  in  order  to  get  a  better  shape  fit 
we  have  to  assume  that  &/d  is  smaller  in  the  N-film  than 
in  the  S-film  so  that  B^/B^  -  10.  This  is  qualitatively 
consistent  with  the  experiment  (see  Chapter  V.3),  but 
because  of  unknown  dependence  of  B  on  £/d  we  cannot  argue 
that  the  ratio  B  /Bc  ^  10  is  reached. 

However  the  experimental  curve  may  be  a  combination 
of  different  McMillan's  curves  such  as  would  be  obtained 
if  the  transmission  probability  a  was  not  constant  over  the 
area  of  the  junction. 

The  temperature  dependence  of  McMillan's  peak 
structure  should  follow  the  curves  marked  by  PI  and  P2  in 
Fig. III. 3.  This  dependence  is  combined  in  the  tunneling 
experiment  with  the  BCS  temperature  dependence  of  the  gap 
edge  in  the  aluminum  electrode,  to  form  the  temperature 
dependence  of  PI  ±  A  and  P2  ±  A.  The  results  of  a  typical 
experiment  are  plotted  in  Fig. VI. 8.  At  most  temperatures 
the  structure  P2  -  A  was  not  resolved  but,  since  there  are 
four  structures  and  only  three  unknowns  PI,  P2,  and  A,  it 
is  possible  to  calculate  the  unknowns  at  each  temperature 
from  three  structures  and  reconstruct  the  position  of  the 
fourth  or  check  the  consistency  of  the  method  where  all 
four  structures  are  resolved.  Also  A,  the  gap  of  the 
aluminum  electrode,  can  be  compared  with  the  temperature 
variation  of  a  BCS  gap  with  Tc  =  1.7°K  as  has  been  done 
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FIG.  VI. 8 


The  temperature  dependence  of  the  peak  positions 
0  experimental  points  corresponding  to  PI  -  A  f  pi  +  A , 
and  P2  +  A  ,  E82,  dAl  =  100  dgn  =  2260  ft. 
o  calculated  positions  of  P2  -  A  using  the  experimental 
points.  Note  the  consistency  of  the  results  at  1.63°K 
where  calculated  and  experimental  points  coincide. 

(H  Pi/  1Z3  P2  positions  calculated  from  experimental 
points.  In  the  lower  graph  the  gap  of  the  aluminum 
electrode  is  shown  from  which  tunneling  into  the  proximity 
structure  was  performed  (O)  as  calculated  from 
experimental  points  in  the  upper  part.  For  comparison 
the  BCS  temperature  variation  of  a  gap  is  included. 


Aai  (meV)  STRUCTURE  POSITION  (meV) 
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in  the  lower  part  of  Fig. VI. 8.  Remarkably  good  agreement  is 
obtained.  The  values  of  P2  from  the  McMillan  model 
calculations  lie  within  the  scatter  of  the  experimental 
values;  whereas  the  values  of  PI  lie  about  0.025  meV  lower 
than  the  experimental  curve  at  all  temperatures. 

Other  experimenters  have  possibly  seen  some  evidence 
of  the  multiple  peaks  in  thicker  films.  Claeson  et  al  (1967) 
observed  structure  with  an  aluminum-lead  sandwich  that  they 
assigned  to  the  superconducting  properties  of  the  contacts. 

It  disappeared  at  lower  temperatures,  possibly  because  the 
smaller  peaks  are  more  difficult  to  resolve  there.  Adkins 
and  Kington  (1969)  observed  "unexplained  dips  in  the  density 
of  states  above  the  gap  and  structure  at  zero  voltage"  on 
the  tunneling  conductance  from  aluminum  into  a  copper-lead 
proximity  sandwich  at  1.25°K.  This  is  just  the  kind  of 
structure  to  be  expected  if  the  aluminum  was  superconducting 
and  their  copper-lead  sandwich  had  a  clean  interface  with  a 
high  a  (see  Fig. VI. 7).  Freake  and  Adkins  (1969)  have 
observed  broad  multiple  peaks  in  normal  metal-proximity 
tunneling  done  at  T  -  0.06°K  and  Hauser  (1966)  reported 
structure  near  A  for  lead  when  tunneling  into  an  aluminum- 
lead  sandwich. 

VI. 2. c  Dependence  of  data  on  a  and  thickness 

In  Table  II  the  r  -values  determined  from  experimental 
data  by  the  three  different  methods  discussed  earlier  are 
compared.  Only  the  peak  position  method  was  used  for  the 
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TABLE  II 


Transmission  probability  a  and  as  a  function 
of  the  oxidation  conditions  for  Al-Sn  proximity 
sandwiches  at  1 . 4 °K. (Tunneling  into  the  aluminum  side 


of  a  proximity  sandwich) . 
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samples  of  Sec. Ill  because  neither  nor  the  position  P2 
could  be  accurately  determined  for  these  weakly-coupled 
(low  a)  sandwiches.  The  a-values  were  determined  using 
eq.III.4  and  putting  =  2  which  may  be  an  assumption  of 
limited  accuracy  because  B  probably  rises  rapidly  above 

this  value  if  the  condition  >>  d.T  is  not  fulfilled. 

N  N 

The  difference  among  the  a-values  determined  by  the 
three  methods  arises  partly  because  bulk  values  of  w  and 

1  for  aluminum  were  used  in  the  self-consistent  calculations 
whereas  it  is  well-known  that  different  values  apply  to  thin 
films  (see  Strongin  et  al ( 1965)  )  .These  bulk  values  mainly 
affect  the  relationship  between  n  and  for  the  McMillan 
model  and  thus  affect  the  reliability  of  the  slope  method 
of  determining  F^  because  of  its  dependence  on  T^.  However 
it  is  clear  that  a  decreases  systematically  as  the  amount 
of  oxidation  at  the  interface  increases,  being  in  about  the 
ratio  80:20:1  for  the  samples  of  Sec. I,  II  and  III 
respectively.  This  range  of  a-values  was  achieved  by 
creating  the  following  conditions  during  the  time  interval 
At  between  the  evaporation  of  the  aluminum  and  tin. 

-7 

Sec.  I  samples:  pressure  ~  2  x  10  Torr ,  At  -  1  mm, 
and  a  shield  around  the  sample  was  held  at  liquid  nitrogen 
temperatures . 

-7 

Sec. II  samples:  pressure  ~  5  x  10  Torr,  At  =  1  min. 

Sec. Ill  samples:  specimen  chamber  isolated  from  pumps, 
pressure  >>  10_7Torr  for  At  -  10  min.  All  pressures  were 
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read  on  a  cold-cathode  ionization  gauge. 

The  process  of  trying  to  control  a  is  very  complex 
and  cannot  be  done  with  precision.  In  a  vacuum  of  3  x  10 
Torr  a  monolayer  strikes  a  clean  surface  in  about  one  sec 
so  one  cannot  hope  for  an  ideally  clean  surface  especially 
on  aluminum,  which  oxidizes  readily.  Adsorbed  layers  are 
partially  penetrated  by  the  next  metal  that  is  evaporated 
so  that  interface  properties  will  depend  on  many  factors 
including  the  specific  metals  involved  (see  Handy  (1962)). 

Finally  it  may  be  noted  from  Table  II  that  the 
measurable  parameters  depend  on  the  thickness  d^  in  accord 
with  the  McMillan  model.  Samples  from  different  sections 
cannot  be  compared  because  the  ^-dependence  masks  the 
thickness  effect.  Samples  prepared  at  the  same  time  provide 
the  most  striking  comparison;  for  example  the  E  samples 
have  dN (E82) /dN (E83)  =1.5  whereas  from  the  peak  separation 
method  (E83 )  /r  (E8 2 )  =  1.55.  As  expected  the  properties 
are  only  slightly  dependent  on  the  thickness  of  the  tin 
film. 


Experiments  with  thicker  films  would  be  interesting 
but  it  will  then  be  more  difficult  to  satisfy  the  condition 
£  >  d  necessary  for  the  McMillan  model.  It  would  be  valuable 
to  have  a  model  that  would  connect  the  De  Gennes  (1964a) 
assumption  of  the  dirty  limit  in  thick  S-films  with 
McMillan's  assumptions  in  thin  N-films. 


- 


•  3  Tunneling  into  the  S-side  of  a  Proximity 
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Sandwich 


VI. 3. a  Tunneling  from  a  normal  metal 

In  contrast  to  tunneling  into  the  N-side  of  a 
proximity  sandwich  tunneling  from  a  normal  metal  into  the 
S-side  of  a  sandwich  does  not  offer  any  information  that 
will  aid  in  the  identification  of  the  parameters  of  the 
proximity  sandwich.  The  dependences  of  [  1  -  a(0) ]  on 
(  1  -  T/T  ) ,  as  calculated  from  the  McMillan  model,  are 
very  close  to  the  BCS  dependences.  These  results  are  to 
be  expected  because  in  the  S-side  of  the  proximity  sandwich 
the  peak  P2  in  the  density  of  states  is  more  dominant.  For 
all  a  the  temperature  dependence  of  this  peak  (see  Fig. III. 4) 
resembles  the  temperature  dependence  of  the  gap  in  a  BCS 
superconductor  and  therefore  the  behaviour  of  the  proximity 
sandwich  must  be  BCS-like. 

The  observations  agree  quite  well  with  the 
theoretical  predictions.  The  plot  of  the  experimental 
dependences  for  three  junctions  with  different  a  is 
presented  in  Fig. VI. 9.  The  slopes  are  very  close  to  unity 
but  in  some  cases  the  lines  are  shifted  toward  smaller 
values  of  1  -  a(0) .  The  largest  shift  was  observed  for 
junction  030  in  Fig. VI. 9.  The  shift  could  be  caused  by 
a  background  conductivity  (zero  bias  anomaly)  that  varies 
v;ith  temperature,  or  by  variation  of  ot  over  the  sample 
area ,  or  by  a  combination  of  both  effects.  A  variation  of  ot 


t  v.--?  mQ*l **  *^SPG|^^H 
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FIG.  VI.  9 

The  experimental  dependence  of  [  1  -  a (V  =  0) ]  on 
(1  -  T/T^) .  a (V  =  0)  is  the  normalized  tunneling  conductance 
at  zero  bias  between  a  normal  metal  and  the  S-side  of  an 

o  1  i  vm  t  m  vm  —  4-  i  n  nrrwi  mi  f  tt 
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could  affect  the  plot  in  the  following  way.  The  measured 
critical  temperature,  ,  will  correspond  to  the  part  of 
the  junction  with  the  smaller  a  and  therefore  1  -  T/T 
will  be  shifted  to  higher  values,  while  cr(0)  will  be  equal 
to  that  for  the  majority  of  the  junction  area. 

VI. 3. b  Tunneling  from  a  superconductor 

The  normalized  conductance  for  tunneling  from  a 
superconductor  into  the  S-side  of  a  weakly  -  coupled 
proximity  sandwich  exhibits  the  multiple  -  peaked  structure 
shown  in  Fig. VI. 10. b.  The  curve  was  constructed  for 
tunneling  at  1.1°K  between  the  density  of  states  for 
aluminum  (see  Chapter  VI. 1)  and  the  McMillan  density  of 
states  shown  in  Fig. VI. 10. a.  The  expected  positions  of  the 
structure  are  indicated  by  arrows  at  PI  ±  A  and  P2  ±  A.  One 
can  see  very  poor  coincidence  between  the  expected  energy  of 
the  structure  and  actual  energy  of  the  peak  maxima.  Also,  in 
certain  cases  a  peak  appears  in  the  normalized  tunneling 
conductance  associated  with  a  dip  in  the  McMillan  density 
of  states  so  that  extra  peaks  sometimes  appear.  These 
effects  make  final  interpretation  of  the  tunneling 
conductance  curves  quite  difficult. 

In  Fig. VI. 11  an  experimental  curve  at  1.4°K  for 

junction  022  is  compared  with  a  calculated  curve  using 

T  =0.03  meV  and  T  =  0.012  meV .  The  value  for  rg  was  found 
S  N 

by  fitting  the  peak  position  to  Fig. III. 7  and  was  found 

from  the  condition  =  stained  by  assuming 


- 
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FIG.  VI.  10 


Illustration  of  the  occurrence  of  multiple  peak 
structure  in  the  tunneling  conductance  between  an  ordinary- 
superconductor  and  the  S-side  of  an  aluminum-tin  proximity 
sandwich . 


(a)  The  normalized  density  of  states  in  the 

superconducting  aluminum,  (characterized  by 

A  =  0.24  meV ,  <5^  =  0.013),  and  proximity 


o^n«^T.7T  ("•n  .  N  ( phs  ra  rfpri  T  =  0.05  rrif^V  . 

- -  --prox  '  S  . 1  ' 

Tn  =  0.01  meV,  A|h  =  0.539  meV ,  and  A?h  =  0.141 

meV) .  When  voltage  is  applied  the  densities  of 

states  shift  on  the  energy  scale, 

(b)  The  resulting  calculated  tunneling  conductance 

at  1 . 1°K  ( - )•  Expected  structures  at  PI  ±  A 

and  P2  ±  A  are  indicated  by  arrows. 

For  comparison  the  tunneling  conductance  between  a 


normal  metal  and  the  same  proximity  sandwich  at  the  same 
temperature  is  also  plotted  (  ) . 


-0.5 


0.5 
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FIG.  VI.  11 


Comparison  of  experimental  and  theoretical 
normalized  conductance  for  tunneling  from  superconducting 
aluminum  into  the  S-side  of  aluminum-tin  proximity  sandwich. 
-  experiment  022,  T  =  1.4°K,  d^  =  2790  8, 


d_  =  1015  8. 

Sn 

- calculated,  F<-,  =  0.03  meV ,  F^  =  0.012  meV , 

A^h  =  0.55  meV,  A^h  =  0.10  meV. 

S  N 

Aluminum  electrode  from  which  the  tunneling  was  done  is 

=  0.12  meV ,  and  6  = 


characterized  by  A 


0.013. 
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BN  “  BS  and  VFN  =  VFS*  The  Peak  positions  agree  quite  well, 
but  the  structure  on  the  calculated  curve  is  sharper  than 
the  experimental  structure.  However  the  experimental  curve 
may  be  a  combination  of  different  McMillan  curves  such  as 
would  be  obtained  if  the  transmission  probability  a  was  not 
constant  over  the  area  of  the  junction. 

The  temperature  dependence  of  PI  ±  A  and  P2  ±  A 
should  arise  from  the  combination  of  the  temperature 
dependence  of  the  McMillan  peak  structure  (see  Fig. III. 4) 
and  the  BCS  temperature  dependence  of  the  gap  edge  in  the 
aluminum  electrode.  The  results  of  a  typical  experiment  are 
plotted  in  Fig. VI. 12.  Since  there  are  four  structures  and 
only  three  unknowns  PI,  P2,  and  A,  it  is  possible  to 
calculate  the  unknowns  at  each  temperature  from  three 
structures  and  reconstruct  the  position  of  the  fourth  or 
check  the  consistency  of  the  method  where  all  four 
structures  are  resolved.  The  values  of  PI  and  P2  as 
calculated  from  the  experimental  points  of  PI  ±  A  and 

P2  ±  A  are  compared  with  the  theoretical  values  ( - ) . 

Also  A,  the  gap  of  the  aluminum  electrode,  can  be  compared 
with  the  temperature  variation  of  a  BCS  gap  with  T  =  1.6°K 
as  has  been  done  in  the  lower  part  of  Fig. VI. 12.  Remarkably 
good  agreement  is  obtained. 

VI. 3. c  Dependence  of  data  on  a 

In  Tables  III  and  IV  the  ^-values  determined  from 
the  experimental  data,  and  the  transmission  probabilities  a 
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FIG.  VI.  12 


The  experimental  temperature  dependence  of  the  peak 
positions . 

^  experimental  points  corresponding  to  A-  PI,  PI  +  A , 
P2  -  A,  and  P2  +  A,  junction  022. 
o  calculated  values  using  the  experimental  points. 
Note  the  consistency  of  the  results  for  A  -  PI. 

El  pi  ,  D  P2  positions  calculated  from  experimental 
points . 

In  the  lower  graph  the  gap  of  the  aluminum  electrode 
is  shown  from  which  tunneling  into  the  proximity  structure 
was  performed  (O)  as  calculated  from  the  experimental  points 
in  the  upper  part.  For  comparison  the  BCS  temperature 
variation  of  a  gap  is  included. 
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calculated  from  them,  are  presented.  The  r  values  were 
determined  by  two  different  methods;  first  by  comparison  of 
the  experimental  critical  temperature  T  with  Fig.  III. 5, 
and  secondly  by  comparison  of  the  experimental  values  of 
peak  positions  with  Fig. III. 7.  The  a~values  were  determined 
using  eq.III.4  and  putting  B  =  2. 

u 

In  Table  III  there  are  presented  the  junctions  in 
which  the  transmission  probability  a  was  controlled  by 
different  amounts  of  oxidation  of  the  tin  layer. xt  is 

clear  that  o t  decreases  systematically  as  the  amount  of 
oxidation  at  the  interface  increases,  being  in  about  the 
ratio  200:30:1  for  the  samples  of  Sec. I,  II,  and  III 
respectively.  This  range  of  a  values  was  achieved  by  creating 
the  following  conditions  during  the  time  interval  At 
between  the  evaporation  of  tin  and  aluminum. 

-7 

Sec. I:  pressure  ~5  x  10  Torr ,  At  -  60  sec. 

Sec. II:  pressure  ~0.2  Torr  of  oxygen  dried  in 
liquid  nitrogen  trap,  At  -  100  sec. 

Sec. Ill:  pressure  ~0.2  Torr  of  dry  oxygen, 

At  -  1000  sec. 

When  Sec. I  of  Table  III  is  compared  with  Sec. II  of  Table  II 
(in  which  the  oxidation  conditions  are  about  the  same)  one 
can  immediately  see  that  the  a  obtained  by  the  oxidation 
of  aluminum  is  much  smaller  than  the  cl  obtained  by  the 
oxidation  of  tin,  which  is  to  be  expected  because  aluminum 

oxidizes  more  readily  then  tin  does. 

The  a-value  derived  from  the  critical  temperature  is 


■ 


TABLE 


III 


Transmission  probability  a  and  F  as  a  function 

O 

of  the  oxidation  conditions  for  Sn-Al  proximity  sandwiche 
Peaks  PI  and  P2  are  measured  at  1.4°K,  unless  otherwise 
stated.  (Tunneling  into  the  tin  side  of  a  proximity 
sandwich) . 
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smaller  than  that  derived  from  the  peak  positions.  The 
difference  may  be  explained  by  slight  non-uniformity  of 
the  oxide  layer  which  forms  the  barrier  between  the  films 
in  proximity.  Regions  with  thicker  oxide  will  have  a 
smaller  value  of  a  and  higher  value  of  T  .  In  the  experiment 
we  will  detect  the  highest  in  the  sample,  thus  getting  a 

smaller  value  of  a  than  the  average  value  which  is  measured 
by  the  peak  positions. 

The  discrepancy  between  the  a  values  determined  by 
the  two  methods  is  most  remarkable  in  Sec. I.  This  probably 
arises  from  the  fact  that  the  McMillan  model  is  valid  only 
for  a  <<  1,  while  the  values  of  a  in  this  section  range 
from  0.15  to  0.70. 

We  have  also  attempted  to  control  a  by  another  means. 
The  results  are  collected  in  Table  IV. 

A  thin  layer  of  SiO  was  evaporated  in  between  the 
aluminum  and  tin  in  junctions  K  and  L.  The  thickness 

of  SiO  was  calculated  from  the  density  p  =  2.165  g/cm3 
and  known  mass  sensitivity  of  the  frequency  of  the  quartz 
crystal  thickness  monitor.  The  sandwiches  produced  in 
this  way  usually  exhibited  more  peaks  in  the  conductance 
curve  than  predicted  by  the  McMillan  model.  When  the  peaks 
were  plotted  against  temperature,  it  looked  as  if  peaks 
from  different  junctions  were  drawn  onto  the  same  sheet 
of  paper.  We  assigned  this  behaviour  to  the  inhomogeneity 
of  the  SiO  layer,  which  probably  condenses  in  islands.  The 
resulting  effect  is  equivalent  to  two  junctions  in  parallel, 
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TABLE  IV 


Transmission  probability  a  for  barriers  formed 
by  the  evaporation  of  SiO  or  glow  discharge  oxidation. 
Peaks  PI  and  P2  are  meassured  at  1.4°K.  (Tunneling  into 
the  tin  side  of  a  proximity  sandwich) . 
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The  glow  discharge  current  was  1mA  for  60  sec  in  0.15  Torr  of  dry 
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one  with  a  high  value  of  a  and  the  other  one  with  a  small 
value  of  cl.  Also  a  very  large  discrepancy  between  a 
determined  from  Tc  and  from  the  peak  positions  (the 
strongest  sets  of  peaks  were  chosen)  speaks  for  a  large 
variation  of  a. 

The  glow  discharge  method  was  employed  to  control 
the  value  of  a  in  the  junctions  I  and  J.  The  method 
produced  very  thick  oxide  and  resulted  in  very  weak  coupling 
between  the  films. 
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CHAPTER  VII 

SUGGESTIONS  FOR  FURTHER  WORK 

A  most  interesting  extension  of  the  present 
experiments  would  be  the  fabrication  of  a  double- junction 
in  which  the  tunneling  characteristics  of  a  proximity 
sandwich  could  be  measured  from  both  aluminum  and  tin 
sides  simultaneously.  One  should  then  be  able  to  compare 
the  value  of  a  and  the  positions  of  the  peaks  PI  and  P2 
as  measured  from  both  sides  of  the  sandwich.  However  such 
a  double- junction ,  requiring  the  evaporation  of  four  metal 
films  and  three  oxidation  procedures,  would  be  rather 
difficult  to  fabricate. 

Of  some  interest  would  be  experiments  in  which  the 
critical  temperature  of  the  base  layer  was  higher  than  the 
critical  temperature  of  the  proximity  sandwich  (e.g. 
junctions  Pb-I-Af,  Sn ,  or  Sn-I-Sn,  A1 ,  where  I  stands  for 
the  tunneling  barrier) .  One  should  be  able  to  measure  the 
temperature  dependence  of  the  peak  positions  over  a  broader 
range  of  reduced  temperature  T/Tc  in  such  a  junction.  Also 
the  peaks  would  be  further  separated  and  they  would  not 
be  mixed  up  at  V  =  0. 

The  McMillan  theory  should  also  be  tested  for  a  pair 
of  metals,  one  of  which  is  a  superconductor  and  the  other 
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one  a  "real"  normal  metal,  which  is  known  not  to  be 
superconducting .  The  best  choice  would  probably  be  copper 
and  lead,  which  have  already  been  extensively  studied  in 
the  De  Gennes  approximation  (i.e.  thick  films  with  short 
mean  free  path  and  with  coupling  between  them  as  strong  as 
possible) . 

From  the  theoretical  point  of  view  the  general 
dependence  of  the  McMillan  model  on  the  mean  free  path 
should  be  established,  or  at  least,  introduction  of  the 
diffusion  limit  into  the  McMillan  model,  as  suggested  by 
Deutscher  and  De  Gennes  (1969)  ,  should  be  done.  This  should 
open  further  possibilities  for  the  experimental 
investigation. 
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CHAPTER  VIII 


CONCLUSION 


The  tunneling  characteristics  of  junctions  between 
aluminum  and  both  sides  of  proximity  sandwiches  formed  by 
evaporated  films  of  aluminum  and  tin  have  been  measured 
and  the  experimental  results  compared  with  a  computer 
analysis  using  the  McMillan  model. 

To  fulfil  the  assumptions  of  the  McMillan  model  the 
coupling  of  the  proximity  films  has  been  limited  by  allowing 
slight  oxidation  at  the  interface.  The  oxidation  of  the 
interface  in  dry  oxygen  or  in  "vacuum"  (pressure  -  10  ^Torr) 
produced  values  of  the  transmission  probability  a  which  are 
qualitatively  quite  consistent  with  the  different  oxidation 
conditions  and  the  kind  of  oxidised  surface  (either  tin  or 
aluminum) .  Glow-discharge-assisted  oxidation  produced 
extremely  weak  coupling  between  the  metals  in  proximity. 

Thin  evaporated  layers  of  SiO  were  not  suitable  for 
decoupling  of  the  metals,  because  they  were  not  homogeneous. 

The  junction  conductance  at  zero  bias  when  tunneling 
from  normal  aluminum  into  the  aluminum  side  of  the  proximity 
sandwich  has  been  compared  in  the  vicinity  of  the  critical 
temperature  of  the  sandwich  with  the  empirical  expression 

1  -  a(0)  =  constant  (1  -  T/T^)n. 


VIII. 1 


' 


n 
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Values  from  n  =  1  to  6  have  been  obtained.  The  McMillan 
model  predicts  the  observed  magnitude  of  a (0)  and  very 
nearly  this  temperature  dependence,  although  for  the  higher 
n-values  the  experimental  data  fit  the  empirical  law 
somewhat  better  than  the  model  calculations  do.  The 
measured  conductance  as  a  function  of  voltage  fits  the 
McMillan  model  very  well. 

When  tunneling  from  normal  aluminum  into  the  tin 
side  of  the  sandwich  is  performed  in  the  vicinity  of  the 
critical  temperature  of  the  sandwich  the  McMillan  model 
predicts  the  BCS  value  n  =  1  in  relationship  VIII. 1 
for  the  junction  conductance  at  zero  bias.  Experimentally 
n  =  1  also  but  sometimes  the  constant  required  to  fit 
eq.VIII.l  is  smaller  than  for  BCS  superconductors.  The 
shift  in  constant  may  be  attributed  to  a  variable  background 
conductivity  and  to  inhomogeneity  of  the  decoupling  barrier. 

When  tunneling  from  a  superconductor  into  a 

proximity  sandwich  is  performed  a  multi-peaked  conductance 
curve  is  observed  and  compares  well  in  some  detail  with  the 
predictions  of  the  McMillan  model.  The  positions  of  the 
peaks  and  their  variation  with  temperature  as  well  as  their 
dependence  on  the  proximity  film  coupling  and  thickness  of 
the  normal  (aluminum)  film  have  all  been  compared  with 
theory  and  agreed  well.  Those  discrepancies  that  exist  may 
be  largely  a  result  of  applying  calculations  using 
properties  of  bulk  aluminum  to  experiments  with  thin  films 
for  tunneling  into  the  aluminum  side  of  the  sandwich  or  a 
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result  of  slight  inhomogeneities  in  the  decoupling  barrier 
for  tunneling  into  the  tin  side  of  a  sandwich. 

The  McMillan  theory  predicts  sharper  peaks  than  the 
observed  ones.  The  broader  peaks  in  the  experiment  may  be 
explained  by  superposition  of  different  McMillan  curves 
with  slightly  different  a. 

As  a  whole  the  McMillan  model  accounts  well  for  the 
observed  behaviour  of  aluminum-tin  proximity  sandwiches. 


' 
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APPENDIX  I 


THE  CALCULATION  OF  THE  FUNCTION  B 


The  value  of  the  function  B  in  the  McMillan 
tunneling  model  of  the  proximity  effect  will  be  calculated 
under  the  following  assumptions : 

(a)  the  mean  free  path  is  infinite/ 

(b)  the  reflection  on  the  back  face  of  the  film  is  specular, 

(c)  the  probability  dw  that  an  electron  penetrates  the 
barrier  to  enter  the  film  within  a  solid  angle  d^  at  an 
angle  0  to  the  normal  obeys  the  cosine  law. 

The  situation  is  illustrated  in  the  figure  below 
where  the  plane  of  the  thin  film  of  the  thickness  d  is 
perpendicular  to  the  z  axis. 


dw  =  w0cos9  dft, 

where  dfi=2Trsin0d0 .  w0  can  be  found  from  the 
equation 


tt/2  tt/2 

1  =  /  dw  =  2ttw0  /  cos0sin0  d0 
0  0 


normalization 


=  ttw  0  . 


•  * 
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The  path  which  is  travelled  in  the  film  of  thickness  d  by  an 
electron  which  enters  the  film  at  angle  ©  is 


2MQ) 


2d 

COS0 


Therefore  the  average  path  L  is 


L 


tt/ 2 
=  / 

0 


COS0 


2d 

COS0 


sin0 


d0 


4d  . 


When,  following  McMillan,  we  write  L=2dB,  we  may  immediately 
see  that  B=2.  However,  the  calculation  is  still  valid,  if 
instead  of  assumption  (b)  we  assume  that  the  electron  is 
reflected  at  the  back  face  of  the  film  according  to  the 
cosine  law.  Therefore  we  may  conclude  that  B=2  as  long  as 

3r-w^  t  i~<  ^  P5  rirl  -HT-io  rpf  1  or"|-i  nn  P  +■ 

UUtJ  L  v— •  -i-  ^  i-  -L  »  \  *—4.  /  \  - - - -  ' — -  — -  —  ■ —  “ — -  — 

the  back  face  of  the  film  is  a  mixture  of  specular  and 


diffuse  reflection. 
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APPENDIX  II 


THE  CALCULATION  OF  THE  CRITICAL  TEMPERATURE  OF  THE 


MCMILLAN  MODEL 


The  method  for  calculation  of  the  critical 
temperature  of  the  proximity  sandwich,  as  suggested  by 
Mr.  R.  Teshima,  is  described.  The  method  utilizes  the 
fact  that  goes  to  zero  when  T  goes  to  T^  .  First,  new 

variables  x,  y,  and  Y  are  defined  by  equations 


nh 


xA-  , 


A. II. 1 


The  new  variables  are  introduced  into  eqs. III.l  and  the 
eauations  are  evaluated  in  the  limit  Ag  ->  0  ,  giving 


x 


y 


(i  + 


) 


-Y. 


A. II. 2 


similar  procedure  is  applied  to  eqs  III. 2,  giving 


0) 


N 


=  A„  /C Re(  £  )  tanh  ^  dE, 


Y 


N 


0 


E 


2T 


A. II. 3 
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1 


S 

coc 

Xg  /  Re(  ^  )  tanh  ~  dE . 


A. II. 3 


Eqs .  A.  1 1 . 2  can  be  solved  for  x  and  y  and  the  results  put 
into  eqs. a. II. 3 , giving 


Y 


1 


w”  yE2+(rN  +  yrc)  (r„  +  re) 


xn  f 
o 


N 


E2  +  (r  +  r0) 2 

n  s 


—  tanh  ~  dE, 


A. II. 4 


10 

Ac  / 
b  0 


c  E^+  <rN  +  Yrs>  <rN  +  FS> 


+  (PN  +  V 


tanh  2^  dE . 


When  y  is  eliminated  from  eqs.A.II.4  we  c3e t  the  final 
result 


f  = 


I  - 


v4 

-  _ _  .  _ -  +  +  Aer..i.  -i  =  o 

ANI2/U-+r^)  "  b  w  J 


XNASrNrsIlI3 


N  S 


where 


N  S  1 


N 


A. II. 5 


N 


uc 

I  =  /  D (E)  dE, 
1  0 

S 

L0C 

I  =  /  D(E)  dE, 
J  0 


N 

(J0C 

I  =  /  E2D(E)  dE, 
Z  0 

S 

u)c 

I  =  /  E2D(E)  dE, 
0 


and 


r  +  Tg  g 

D(E)  =  e2  +  (r  +  rc)2“  e  tanh  2t  * 


N  '  *  s 

The  temperature  for  which  F  =  0  is  the  critical 


temperature  of  the  McMillan  model.  The  equation  F  =  0  has 
more  than  one  solution.  The  root  of  physical  interest  is 
selected  from  the  temperature  interval  1.2°K  to  3.8°K. 


' 
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GLOW  DISCHARGE  OXIDATION  OF  FRESHLY  EVAPORATED  ALUMINUM 

FILMS 


Introduction 

A  discussion  is  presented  of  glow  discharge  assisted 
oxide  growth  on  the  surface  of  freshly  evaporated  aluminum 
films  at  14°C  with  special  attention  to  the  preparation  of 
tunnel  junctions.  The  mass  gain  of  the  aluminum  film  during 
the  oxidation  was  measured  with  a  quartz  crystal  microbalance 
and  recorded  as  a  function  of  time. 


m  /<>  n  4-  /-n  r-v  -ft 

nit:  x_  ex  e-vx  w  j_ 


■***''  ^  ~  v-v  4-  d  n  o  nv  -I  Lv  ^  ,^3  V-\  t  r 

iituoo  y  axii  v^aniiu  u  ucov^x  xucu  iJ y 


CX  H  Xii^j  -1.V— 


oxidation  law  covering  the  whole  period  of  oxidation.  The 
oxidation  rate  was  found  to  depend  only  on  the  total  current 
through  the  glow  discharge;  the  applied  voltage  and  oxygen 
pressure  were  not  important. 


Experimental 

The  vacuum  apparatus  was  a  standard  system  used  for 
preparation  of  thin  film  tunnel  junctions  (three  stage  oil 
diffusion  pump  with  diameter  of  4  inches  with  a  liquid 
nitrogen  trap  backed  by  a  mechanical  pump) .  Before  the 
evaporation,  the  vacuum  reached  1  x  10  7  Torr  and  was 
kept  in  the  range  of  10~6  Torr  during  the  evaporation. 

The  quartz  crystal  microbalance  consisted  of  two  5  MHz 
quartz  crystals  mounted  in  the  same  water- cooled  copper  holder. 
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One  crystal  was  directly  exposed  to  evaporated  material  and 
the  surrounding  oxidizing  atmosphere,  the  other  was  in 
thermal  contact  with  the  copper  holder  only.  The  change  of 
beat  frequency,  which  arose  from  mass  deposition  onto  the 
exposed  crystal,  was  recorded.  The  mass  sensitivity  was 
0.02  yg/cm2  Hz.  The  drift  of  the  system  was  about  2  Hz 
per  hour  after  an  initial  warm-up  period. 

A  film  of  aluminum  about  300  A  thick  (99.999%  purity) 
was  evaporated  onto  the  surface  of  the  quartz  crystal  from  a 
tungsten  coil  at  a  distance  of  about  45  cm.  A  few  minutes 
after  the  deposition  was  completed,  commercial  grade 
oxygen  purified  by  passing  through  a  coil  immersed  in  liquid 
nitrogen,  was  admitted  to  the  system  at  various  pressures 
ranging  from  0.02  to  0.2  Torr. 

A  glow  discharge  was  initiated  in  accordance  with  Miles 
and  Smith  by  applying  a  negative  d.c.  potential  to  an  aluminum 
electrode.  The  metallic  part  of  the  evaporator  served  as  the 
ground  electrode.  Current  I  flowing  through  the  negative 
glow  electrode  was  measured  and  its  value  controlled  either 
by  varying  the  applied  potential  or  oxygen  pressure  in  the 
evaporator . 

The  negative  glow  electrode  consisted  of  two  straight 
pieces  of  aluminum  wire  about  3  in  long  and  5  in  apart.  They 
wqtq  connected  together  and  brought  to  the  vacuum  feed-through 
by  a  wire  surrounded  with  10  mm  glass  tubing  in  order  to 
restrict  the  discharge  to  the  aluminum  wires  only.  The 
position  and  the  shape  of  the  glow  electrode  were  chosen 
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to  produce  an  oxide  film  with  uniform  thickness  over  an 

2 

area  of  1  x  3  in  .  This  was  the  size  of  microscope  slide 
on  which  three  tunneling  junctions  were  prepared  at  the 
same  time.  It  was  assumed  that  uniformity  of  the  tunneling 
resistance  of  the  oxide  implies  uniformity  of  the  oxide 
thickness.  First,  it  was  found  experimentally  that  the 
tunneling  resistance  of  an  oxide  film  formed  in  a  glow 
discharge  depends  approximately  on  r  where  r  is  the 

radial  distance  from  a  single  straight  electrode.  Secondly, 
two  such  electrodes  were  combined  to  produce  uniform  oxide 
over  the  desired  area.  The  quartz  crystal  microbalance  was 
placed  in  an  asymmetrical  position  about  12  cm  above  one 
electrode.  The  actual  oxidation  rate  depends  on  the  exact 
geometry  of  the  ground  electrode  and  it  was  found  to  be 
different  in  different  evaporation  systems. 

Tiny  sparks  that  sometimes  appear  on  the  glow  electrode 
may  produce  pinholes  in  the  oxide  and  certainly  introduce 
irreproducibili ty  in  tunnel  junction  preparation.  It  is 
therefore  necessary  to  clean  the  electrode  prior  to  the  actual 
oxidation.  This  may  be  accomplished  when  the  system  is  being 
evacuated  for  the  first  time  by  maintaining  a  glow  discharge 
with  a  higher  voltage  than  planned  for  the  actual  oxidation, 
until  all  sparks  disappear.  In  this  process,  all  dust 
particles  and  surface  irregularities  which  may  form  sharp 
points  on  the  surface  of  the  electrode  are  burned  off. 
Sometimes,  when  aluminum  wires  of  low  purity  are  used,  the 
sparking  cannot  be  avoided. 
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Results  and  discussion 

Typical  oxide  growth  characteristics  are  shown  in  Fig. 
A.III.l.  The  change  of  crystal  frequency  (counting  was 
started  after  the  evaporation  of  the  aluminum  film  was 
completed)  is  plotted  vs  time  for  glow  discharge  oxidation 

at  an  applied  voltage  of  900  volts  and  current  of  3.5  mA 

through  the  glow  electrode.  There  are  three  stages  of  glow 
discharge  oxidation  visible  in  Fig. A.III.l. 

1)  When  the  system  was  filled  with  dry  oxygen  (at  time  ti) 
the  frequency  of  the  crystal  increased  about  Af  (Af  is 
usually  from  4  to  6  Hz) .  This  process  was  usually 
completed  within  50  sec.  This  change  did  not  disappear 
when  the  system  was  again  pumped  to  high  vacuum  but  it 

was  then  stable  in  time.  It  was  concluded  that  the 

frequency  jump  Af  corresponds  to  a  formation  of  aluminum 
oxide  with  its  limiting  thickness,  as  discussed  by  Mott 
(1939,1940,1947),  Mott  and  Cabrera  (1949),  Fromhold  and 
Cook  (1967),  and  Hart  (1956).  To  convert  mass  gain  data 
per  unit  area  to  change  in  oxide  film  thickness,  we  note 
according  to  Dignam  et  al.  (1966)  that  a  mass  gain  of 
1  yg/cm2  corresponds  to  an  increase  of  oxide  thickness  of 
66.4  8.  Therefore  a  frequency  change  of  5  Hz  corresponds 
to  0.1  yg/cm2  or  7  8  of  aluminum  oxide.  This  seems  to  be 
a  bit  too  low  for  the  saturated  thickness  of  aluminum 
oxide  (which  is  usually  about  25  8)  but  oxide  formation 
on  the  surface  of  very  thin  evaporated  films  in  dry 
oxygen  is  probably  different  from  that  on  bulk  aluminum. 
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FIG.  A. Ill . 1 


The  change  of  the  crystal  frequency  (  left  hand  scale  ) 
or  corresponding  mass  gain  (  right  hand  scale  )  plotted 
vs  time  (  glow  discharge  oxidation  at  I  =  3.5  mA  )• 

Q  experimental  points, 

-  least  squares  fit  to  generalized  Mott-Cabrera  law, 


linear  law  fit  • 
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2)  When  the  glow  discharge  was  switched  on  (at  the  time  tz) 
the  frequency  steadily  increased  with  time  (which 
represents  the  forced  growth  of  the  aluminum  oxide) . 

This  stage  is  discussed  later  on  in  more  detail. 

3)  When  the  glow  discharge  was  switched  off  (at  the  time  t3) 
the  frequency  usually  dropped  about  2  Hz.  This  stage 
probably  corresponds  to  the  outgassing  of  the  loosely 
bound  or  physically  adsorbed  oxygen  molecules  from  the 
surface.  When  the  system  was  again  pumped  to  high  vacuum 
the  frequency  of  the  crystal  microbalance  remained 
unchanged . 

The  oxide  thickness  that  produced  an  Al-Al2C>3-Sn  tunnel 
junction  with  a  resistance  of  100ft  was  estimated  to  be  about 
27  8  or  about  0.4  yg/cm2  of  mass  gain.  Using  the  known  time 
for  the  junction  preparation,  the  geometrical  factor  was 
employed  to  calculate  the  corrected  time  at  the  position  of 
the  crystal  and  the  oxide  thickness  was  read  from  the 
corresponding  mass  gain  vs  time  plot.  One  has  to  remember 
that  the  oxide  thickness  of  27  8  is  reduced  as  the  top 
electrode  is  evaporated  and  the  metal  ions  penetrate  the  oxide 
As  the  aluminum  oxide  grows , aluminum  ions  must  traverse 
the  oxide  film  to  the  oxide-oxygen  interface  (ionic  current) 
and  simultaneously  electrons  must  be  transferred  through  the 
oxide  (electronic  current)  to  account  for  the  change  of  the 
charge.  It  is  assumed,  that  in  the  present  case,  the  ionic 
current  is  rate  limiting  and  that  we  can  write  for  the  rate 
of  growth  according  to  Cabrera  and  Mott  (19  49)  ,  Fromhold  and 


' 


Cook  (1967) ,  and  Dignam  et  al.  (1966) 
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where  X  is  the  oxide  film  thickness  at  time  t,  ft  is  the 
volume  of  oxide  per  one  ion,  n  is  the  concentration  of 
mobile  ions,  v  is  the  ionic  attempt  frequency,  a  is  the  ionic 
jump  distance,  w  is  either  the  metal-oxide  interface  barrier 
or  barrier  inside  the  oxide  film  (whichever  is  bigger) ,  y 
is  the  ionic  jump  distance  times  effective  charge  and  E  is 
the  electric  field  in  the  oxide.  We  may  derive  two  limiting 
laws  of  growth. 

1)  When  E  =  C,  where  C  is  constant,  we  obtain  a  linear  law 
of  growth 
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A.III.l 


2) 


When  the  potential  V  = 
we  may  write  according 
Dignam  et  al.  (1966) 

~  =  2A  sinh  (|)  ,  B 


EX  across  the  film  is  constant, 
to  Fromhold  and  Cook  (1967)  ,  and 


yEX  =  yV 
kT  kT 


and 
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A. Ill . 2 


where  X  is  the  oxide  thickness  at  time  t  =  0.  Eq.A.III 
o 

2  is  called  the  generalized  Mott-Cabrera  law. 

In  general,  neither  E  nor  V  are  constants.  For  thermal 
oxidation ,  V  is  normally  considered  to  be  the  initial 
difference  between  the  metal  Fermi  level  and  the  oxygen  0 
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level  (see  e.g.  Fromhold  and  Cook  (1967)) .  However,  in  our 
case,  V  also  includes  the  effect  of  the  glow  discharge. 

A  linear  plot  of  mass  gain  vs  time  (circles)  is  shown  in 
Fig.  A.III.l  for  typical  oxidation  at  constant  current, 

1=3.5  mA .  The  solid  line  represents  a  least  squares  fit 
of  the  data  to  Eq .  A. III. 2  (fitted  by  IBM  computer)  .The  fit  is 
quite  good  for  t  -  t2  >  300  sec  (or  for  thickness  greater 
than  27  8) .  For  t  -  t2  <  300  a  linear  law  fits  the  data 
best.  In  terms  of  eqs.A.III.1  and  A. III. 2  it  appears  that 
in  the  initial  stage  of  a  glow  discharge  oxidation  performed 
at  constant  current  there  is  a  rather  constant  electric  field 
in  the  oxide.  That  is,  as  the  oxide  thickness  increases,  the 
potential  drop  across  it  must  increase.  The  later  oxidation 
stage,  t  -  t2  >  300  sec,  which  fits  the  generalized  Mott- 
Cabrera  law  quite  well  is  characterized  by  a  constant  potential 
V  across  the  oxide  film. 

From  the  least  squares  fit  of  data  to  eq .  A. III. 2  we 

have  found  for  the  later  oxidation  stage  at  14°C  and 

1=1  mA  the  values,  B  =  29.2  8  and  A  =  3.13  x  10  3  8  sec  1 , 

and  for  I  =  3.5  mA,  B  =  68.1  8  and  A  =  1.80  x  10_3  8  sec  l. 

These  values  of  B  lead  to  yV  =  0.72  e8v  for  1=1  mA  and 

yV  =  1.69  eAV  for  I  =  3.5  mA.  Because  y  is  connected  with  the 

oxide  structure  and  should  be  about  the  same  in  both  cases  we 

see  that  the  total  current  I  through  the  glow  discharge  affects 

the  potential  across  the  growing  oxide  film.  We  may  compare 

our  results  with  the  high  temperature  oxidation  data  of  Dignam 

o 

et  al.  (1966)  extrapolated  to  14°C,  which  yield  yV  =  8.2  eAV 
and  A  =  2.2  x  10“22  8  sec"1.  The  differences  indicate  a 
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substantial  difference  of  oxide  structure  for  glow  discharge 
oxidation  and  thermal  oxidation.  There  is  of  course  some 
question  about  whether  conclusions  based  on  high  temperature 
oxidation  can  be  extrapolated  to  low  temperature  and  whether 
the  rate  limiting  process  (barriers  inside  oxide  or  in  metal- 
oxide  interface)  remains  the  same. 

In  Fig.  A. III. 2  the  oxidation  data  for  different  runs  are 
plotted  so  that  a  parabolic  relationship  x2  «  (t  -  t2)  would 
give  a  straight  line.  It  may  be  seen  that  a  parabolic  law 
(dashed  line)  fits  the  data  over  a  more  limited  range  than 
does  Eq.A.III.2  (solid  line). 

In  Fig.  A. III. 3  the  data  are  plotted  in  accord  with 
the  direct  logarithmic  law  x  °c  log  (t  -  tz)  •  The  fit  is 
quite  good  for  t  -  tz  <  1700  sec  (thickness  <  43  $) ,  but  for 
longer  oxidation  times  experimental  data  deviate  from  logarithmic 
law.  The  logarithmical  law  fits  the  data  in  the  range  suitable 
for  the  preparation  of  good  tunnel  junctions,  therefore  one 
can  expect  approximately  a  linear  relationship  between  the 
junction  resistivity  and  oxidation  time. 

A  very  important  feature  of  glow  discharge  oxidation  is 
illustrated  in  Fig. A. III. 2.  The  oxide  thickness  at  any 

time  depends  only  on  the  total  current  I  through  the  glow 
discharge  and  is  independent  of  the  way  it  is  controlled.  The 
glow  discharge  produces  the  same  thickness  of  aluminum  oxide 
at  a  low  pressure  of  oxygen  and  high  applied  voltage  as  at 
high  pressure  of  oxygen  and  low  applied  voltage  as  long  as  the 


- 
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FIG.  A. III. 2 

The  change  of  the  crystal  frequency  vs  time  for 
different  oxidation  runs  plotted  in  terms  of  parabolic  law. 
$  represents  the  oxidation  when  1=3.5  mA,  constant 
applied  voltage  is  900  V  and  I  is  kept  constant  by 
varying  pressure  of  dry  O2  from  0.1  Torr  to  0.2  Torr; 
Q)  represents  the  oxidation  when  I  =  3.25  mA,  constant 
pressure  of  dry  02  is  0.11  Torr  and  I  is  kept  constant 
by  varying  voltage  from  800  V  to  1450  V.  Similarly 
represents  the  oxidation  when  1=1  mA,  constant 
applied  voltage  is  900  V  and  I  kept  constant  by 
changing  the  pressure  from  0.03  Torr  to  0.1  Torr,  and 
n  is  oxidation  when  I  =  mA,  constant  pressure  of  dry  02 
is  0.08  Torr  and  I  kept  constant  by  varying  voltage 


from  600  V  to  800  V. 
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FIG.  A. III. 3 


The  crystal  frequency  vs  logarithm  of  time  for 


the  oxidation  run  with  I 


3.5.  mA . 
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current  I  is  the  same.  No  maximum  in  the  oxidation  rate  at 
a  pressure  of  0.05  Torr  as  reported  by  Miles  and  Smith  (1963a) 
was  found. 

Conclusion 

Oxidation  of  thin  evaporated  aluminum  films  has  been 
studied.  The  time  dependence  of  the  oxide  growth  does  not 
conform  to  a  generalized  Mott-Cabrera  law,  a  parabolical 
law,  or  a  logarithmic  law  over  a  broad  time  range.  We  may 
conclude  from  fitting  the  generalized  Mott  -Cabrera  law,  that 
in  the  first  stage  of  glow  discharge  oxidation  the  potential 
across  the  oxide  film  is  growing,  keeping  the  electrical 
field  in  the  oxide  approximately  constant  and  giving  a  linear 
law  of  oxidation,  while  in  the  advanced  stage  of  oxidation, 
the  potential  across  the  film  is  constant,  leading  to  the 
Mott-Cabrera  oxidation  law. 

The  oxidation  rate  was  found  to  depend  only  on  the  total 
current  through  the  glow  discharge. 
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APPENDIX  IV 

THE  EVALUATION  OF  THICKNESS  FROM  THE  INTERFERENCE  FRINGES 

The  interference  fringes  are  photographed  and  their 
spacings  on  the  film  are  measured  using  an  optical 
densitometer.  The  optical  system  of  the  densitometer  is 
focussed  into  a  sharp  spot  (  the  size  is  smaller  than  the 
width  of  a  fringe), over  which  the  intensity  of  the  light 
passing  through  the  film  may  be  recorded.  The  photograph  of 
the  fringes  is  placed  on  a  two-dimensionally  movable  table 
in  the  focal  plane  of  the  optical  system,  and  the  steps 
in  the  fringes  that  correspond  to  the  edge  of  the 
thickness  change  are  adjusted  so  that  a  line  through  them 
would  be  as  nearly  as  possible  parallel  to  one  direction 
of  the  table  motion.  The  positions  of  the  fringes  on  both 
sides  of  the  steps  are  recorded  and  the  thickness 
calculated  from  them. 

A  reference  frame  may  be  defined  by  the  steps  and 
by  one  of  the  fringes.  However,  the  fringes  are  never 
perpendicular  to  the  line  of  the  steps  and  the  photograph 
is  never  placed  on  the  table  of  the  optical  densitometer 
with  the  fringes  exactly  perpendicular  to  one  direction 
of  the  table  motion.  Thus  an  error  that  depends  linearly 
on  the  measuring  distance  from  the  steps  is  introduced. 
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A  method  which  eliminates  this  error  was  designed  and  is 
described  below. 

To  eliminate  the  error,  we  draw  onto  the  photograph 

equally  spaced  lines  (  called  simply  "the  lines"  )  which 

are  approximately  parallel  to  the  steps.  The  line  which 

is  closest  to  the  steps  is  called  the  X-axis.  Assume  for 

a  moment,  that  X  can  be  placed  exactly  parallel  to  the 

direction  of  the  motion  of  the  table  of  the  optical 

densitometer.  See  Fig . A. IV . 1 . a.  The  perpendicular 

distance  between  the  lines  is  a,  the  angle  between  the 

fringes  and  X-axis  is  a,  the  perpendicular  distance 

between  two  fringes  is  &,  the  perpendicular  height  of  the 

step  at  the  edge  of  the  thickness  change  is  A,  and  A' 

and  £*  are  the  same  as  A  and  &  but  the  distances  are 

measured  in  the  direction  parallel  to  X.  It  is  obvious, 

that  j.  =  A |  t  therefore  the  thickness  d  =(X/2)  (A/£) 

does  not  depend  on  the  direction  along  which  we  measure 

the  distances.  The  numbers  A  ,  A  ,  . . .  A  denote  the 

m  n2  nN 

X  coordinates  of  the  intercepts  of  the  fringes  with  the 
line  which  intercepts  the  Y  axis  at  the  point  na,  and  N 
is  the  total  number  of  fringes.  The  intercepts  of  the 
fringes  with  the  lines  which  cross  the  negative  Y  axis  are 
labeled  by  the  symbol  B.  Therefore  from  measurement  we  get 
the  following  two  sets  of  data 

bi  =  AH'  A12'  Al3'  ***  AlN 
A2  =  A21,  A22,  A  2  3  ,  •  •  •  A2N 
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FIG.  A.IV.l 


(  a  )  The  geometry  of  interference  fringes  , 
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where  K  is  the  number  of  lines  which  are  drawn  in  the 
Y  >0  space  and  L  is  the  same  in  Y<0  space.  (  In  Fig. 
A.IV.l.a  K  =  L  =  3  ). 

Now  the  averaging  procedure  for  A'  and  £'  in  terms 
cf  A*s  and  B's  will  be  described.  First  we  define  the 
vectors 


AAk - 

~k+l 

-  Ak  ' 

A . IV . 2 

AB  J  = 

?2  - 

?i+i  ' 

where  k  =  1,2, 

.  .  .  K 

-  1  and  I  =  1,2,  ... 

L 

-  1. 

When 

the  vectors  A, 

~k 

and 

Bj  are  averaged  over 

k 

and 

2  (it 

means  averaging  each  component  over  the  appropriate  fringe) 
we  get  new  vectors 


AA 


1 

K-l 


K-l 

Z 

k=l 


1 

L-l 


L-l 

Z 

2=1 


A. IV.  3 


AB 
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In  the  ideal  case  (  no  scatter  in  the  data  set  A.IV.l  ) 
all  components  of  vectors  AA  and  AB  are  equal.  If  such 
a  component  is  denoted  AA ,  we  can  calculate  the  angle  a 
from  the  expression 


tga 


a 

AA  ‘ 


We  may  define  the  vectors 

Aj  =  Bj  -  Aj  +  j  (AA  +  AB) , 


A. IV. 4 


A . IV  .  5 


where  j  =  1,  2,  ...  M,  and  M  =  min  (  K,  L  ) .  When  the 

I 

vectors  A j  are  averaged  over  j ,  we  get 


A '  = 


1 

M 


M  , 

Z  A  .  . 


A. IV. 6 


In  the  ideal  case  (  no  scatter  in  the  data  set  A.IV.l  ) 
all  components  of  vector  A'  are  equal  to  A ' . 

To  calculate  £ ' ,  we  form  a  vectors  £j^  and  £®  , 

which  components  are  defined  by 


£ 


A 

k ,  n 


Ak , n+1 


A, 

k  ,n  , 


A . IV .  7 


ll,m  BI ,m+l~  BI,m  , 

where  m,  n  =  1,  2,  ...  N  -  1.  Averaging  over  k  and  2,  and 

both  together  we  get  the  vector 


V  = 


K-l 


k=l 
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K-l 
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Again,  in  the  ideal  case  each  component  of  is  equal  to  i '  . 
.  A ' 

Now  we  form  the  vector  by  division  of  each  component  of 

A'  by  corresponding  component  Z 1  .  Averaging  over  the 
components  of  the  vector  A'/Ji'we  get  the  final  answer 


1 

N-l 


N-l 

E 

i  =  l 


A . IV . 9 


The  above  described  method  for  calculation  of  was 

programmed  together  with  its  standard  deviation  for  the 
IBM/360  time  sharing  system  using  APL/360  language 
(  Falkoff  and  Iverson  (1968)). 

Now  it  will  be  shown,  that  when  the  photograph  is 
placed  with  X-axis  at  an  angle  e  to  the  direction  of  the 
motion  of  the  table  of  optical  densitometer,  the  results 
are  not  affected.  For  simplicity  the  calculation  is  done 
only  for  one  line  drawn  in  each  half  plane  Y>0  and  Y<0 
and  for  one  fringe.  The  situation  is  shown  in  Fig . A. IV . 1 . b . 
The  quantities  which  are  read  when  the  angle  e  is 
introduced  are  denoted  with  a  bar.  By  elementary 
geometrical  considerations,  the  following  expressions 
may  be  derived 


*  -  (Bj.  -  Ax  +  2AA)  ( 1  +  fffr), 

A.  IV.  10 

x  -  <s2  -  V  (1  +  -  <V  V  (1  +  & ' 


where  a,=  a 


e  and  the  meaning  of  the  other  symbols  is 


‘ 
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obvious  from  the  Fig . A. IV .1 . b .  One  can  see  from  eq.A.IV.10 

that  the  ratio  — A i s  independent  of  angle  £  . 

I 

When  the  lines  are  not  equally  spaced,  the  resulting 
error  5d  in  the  calculated  thickness  d  is  of  the  order  of 


6d 


2730 


5a 

£tga 


t 


where  5a  is  the  variation  in  the  spacing  of  the  lines,  £ 
is  usually  about  10  mm  and  ot  is  within  10°  around  90° 

(  which  is  easy  to  adjust  ).  If  e.g.  the  maximum  allowed 
error  5d  is  about  5$,  we  may  estimate  5a  -  0.5  mm.  Such 

an  error  is  very  unlikely  because  the  parallel  lines  are 
separated  about  5  mm  and  are  drawn  through  a  rigid  mask 
which  is  machined  from  1/16  inch  brass  plate  . 
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APPENDIX  V 

THE  CIRCUIT  DIAGRAM  OF  THE  QUARTZ  CRYSTAL  MONITOR 

The  5  MHz  quartz  crystals  are  driven  by  two  identical 
transistor  oscillators,  shown  in  Fig.A.V.l.  The  transistors 
T3  and  T6  act  as  impedance  transformers.  Both  oscillators 
are  mounted  inside  the  high  vacuum  jacket  and  the  high 
frequency  from  both  is  brought  out  on  one  wire  together 
with  the  d.c.  power. 

A  circuit  diagram  of  the  mixer  and  filter  is  shown 
in  Fig.A.V.2.a.  Transistor  T7  acts  as  a  mixer  and  detector. 
Transistor  T8  acts  as  a  buffer  separating  the  active  filter 
Al  from  the  mixing  circuit.  The  filter  does  not  pass 
frequencies  higher  than  250  kHz. 

The  final  amplifier  A2  is  shown  in  Fig.A.V.2.b. 

The  output  beat  frequency  is  connected  to  the  input  of 
a  digital  frequency  counter  (  5214L  Preset  Counter, 

Hewlett  Packard  ) . 
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FIG.  A.V.l 


Oscillators  for  quartz  crystal  monitor  which 
are  mounted  inside  vacuum  space.  All  transistors  are 
International  Rectifier,  TR  20-C. 


I - 1 
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FIG.  A.V.2 

(a)  Mixing,  detecting  and  filtering  circuit  for 
extraction  of  beat  frequency  from  two  quartz 
crystal  oscillators. 

(b)  Output  amplifier  for  frequency  counter. 

All  transistors  are  of  the  tvpe  International  Rectifier 
TR  20  —  C;  operational  amplifiers  are  NS  LM  301  A. 
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APPENDIX  VI 

A  SIMPLE ,  CONTINUOUS  READING  LIQUID  LEVEL  DETECTOR 

A  continuous  reading  liquid  level  detector,  employing 
the  difference  between  the  dielectric  constant  of  gas  and 
liquid,  is  described.  The  construction  and  alignment  are 
very  simple  and  the  liquid  level  is  read  directly  from  a 
panel  meter.  Other  authors  (  Williams  and  Maxwell  (1954), 
Meiboom  and  O'Brien  (1963),  Seguin  and  Leonard  (1966), 
Preshkov,  Vetchinkin,  and  Elkonin  (1969),  Dash  and  Boorse 
(1951),  Jones  (1969),  and  Capocci  and  Gugan  (1967))  have 
described  level  detectors  working  on  the  same  principle, 
but  requiring  expensive  instrumentation  or  sophisticated 
electronics  and  more  complicated  methods  for  level  reading. 

In  the  present  system,  the  change  of  the  dielectric 
constant  is  detected  by  a  coaxial  capacitance  probe  which 
forms  one  arm  of  a  capacitance  bridge  that  is  fed  by  a  low 
frequency  voltage.  The  unbalance  voltage  that  arises  from 
the  liquid  level  change  is  amplified,  rectified  and  read 
directly  on  a  panel  meter.  The  bridge  design  does  not  employ 
carefully  balanced  arrangements  with  guards  because  zero 
adjustment  is  provided  with  a  d.c.  voltage  in  the  last 
amplification  stage.  The  system  may  be  used  for  liquid 
helium,  for  liquid  nitrogen,  or  for  other  liquids  by  simply 
changing  the  amplification. 


. 
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The  level  detecting  probe  consists  of  a  coaxial 

capacitance  CQ  of  length  ZQ  =  20.5  cm  (  along  which  the  level 

&  is  changing  )  and  of  the  capacitance  Cg  formed  by  the 

connections  to  the  measuring  device.  The  total  capacitance 

of  the  probe  is  thus  =  C  +  C  .  The  capacitor  C  is  formed 

by  two  stainless  steel  tubes  of  1/4  in-diam  and  3/16  in-diam 

separated  by  nylon  washers  (see  Fig. A. VI. 1)  The  inner  tube 

is  extended  above  the  liquid  region  by  a  s.s.  tube  of  1  mm 

diam  in  order  to  reduce  capacitance  C  .  Small  holes  are 

s 

drilled  along  the  length  of  the  probe  to  permit  easy  liquid 
access  and  prevent  pressure  oscillations  in  the  probe.  Rigid 
construction  is  important. 

When  length  z  of  the  probe  is  immersed  into  the  liquid, 
the  capacitance  increases  about  AC  =  Cq(&/£o)  (Ae/e  ),  where 
£q  is  the  dielectric  constant  of  the  vapor  and  As  is  its 
difference  from  the  liquid  phase.  (  For  liquid  helium 
Ae/eo  =  0.05,  and  for  liquid  nitrogen  As  /e  q  =  0.45  ). 

The  probe  together  with  capacitors  C^,  C^r  and  form 
a  capacitance  bridge  (  see  Fig.  A. VI. 2  ).  The  bridge 
capacitors  are  zero  temperature  coefficient  units  (  type  TCZ , 
Centralab  )  with  values  C1  =  =  82  pF  and  C2  =  62pF  in 

parallel  with  a  20pF  air  trimmer.  The  bridge  is  fed  from  an 
oscillator  with  a  frequency  of  about  20kHz  and  a  peak  output 
voltage  Eq  =  8V.  Its  stability  and  actual  frequency  are  not 
important.  A  simple  double  integrator  circuit  with  a 
regenerative  feedback  oscillator  (  Burr  —  Brown  (19  63) )  is 


used. 


■ 
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FIG.  A.  VI .  1 


The  level  detecting  capacitance  probe. 


BNC  CONNECTOR 
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FIG.  A. VI. 2 


Circuit  diagram  (  The  capacitance  bridge  and 
differential  amplifier  A1  are  connected  to  the  remaining 
system  by  a  long  cable.  ) 


osc. 
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We  can  write  (  Harris  (1957) )  for  the  unbalance  voltage 
U  of  the  bridge  working  into  infinite  impedance, 

C, 


U  =  E 


a  Z 


o  C-.+C.  1+C  /C  +C ./C  +a£ 


,  where  a  = 


As 


e  Z 
o  o 


For  liquid  helium  when  aZ  <<  D  =  1  +  C  /C  +  CA/C  ,  the 

unbalance  voltage  U  is  linearly  proportional  to  the  liquid 

level  height  Z .  The  magnitude  of  the  voltage  U  is  roughly 

proportional  to  As  (  1  +  Cs/Cq  ) -1  (  because  -  C 3  ).  In 

order  to  keep  U  large  (  As  is  small  for  liquid  helium  )  the 

bridge  circuit  is  mounted  directly  on  top  of  the  level 

detecting  probe  and  the  capacitance  C  is  kept  small  by  the 

s 

probe  construction.  In  our  case  C  /C  -  0.8.  This 
^  s  o 

arrangement  requires  less  overall  amplification  and  therefore 

the  requirement  for  component  stability  is  milder. 

When  the  same  system  is  used  with  liquid  nitrogen,  the 

condition  olZ  <<  D  is  not  as  well  satisfied  and  the  voltage  U 

is  a  slightly  non-linear  function  of  the  level  height.  In 

another  system,  which  was  used  for  liquid  nitrogen  only,  the 

ratio  C  /C  =4,  and  U  was  then  a  linear  function  of  Z  . 
s  o 

The  unbalance  a.c.  voltage  U  is  detected  by  the 
differential  amplifier  (  Burr-Brown  (1963))  Al  (  see  Fig. 

A. VI. 2  )  which  has  an  input  impedance  at  least  ten  times  the 
impedance  of  the  capacitors  in  the  bridge  at  the  operating 
frequency.  The  signal  is  amplified  ten  times  and  rectified 
(  Burr  -  Brown  (1963))  by  A2  and  filtered.  A  small  negative 
voltage,  which  serves  as  a  zero  adjustment  is  added  to  the 
d.c.  voltage  at  the  input  to  the  last  amplification  stage  A3. 


- 
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Therefore  it  is  sufficient  to  bring  the  bridge  into  an 
approximate  balance  and  then  to  zero  the  system  by  the 
potentiometer  PI  (350q)  when  there  is  no  liquid  present. 
Since  the  detection  is  not  phase  sensitive,  one  must  insure 
that  the  bridge  does  not  pass  through  the  exact  balance 
condition  as  the  level  of  liquid  gas  is  increased.  The 
liquid  level  is  read  on  panel  meter  (  1  ma  F . S .  )  and  the 
maximum  deflection  (  when  the  system  is  full  of  liquid  )  is 
adjusted  by  R(lkq)  and  P2  (350i3).  The  values  of  these  two 
components  must  be  changed  if  C  and  C  are  changed.  The 
switch  S  changes  the  amplification  factor  by  ten  and  is  for 
liquid  helium/nitrogen  in  the  position  1/2.  All  operational 
amplifiers  used  are  of  the  type  NS  LM  301  A. 

The  electronics  system  follows  the  bridge  unbalance 
voltage  U  linearly  and  disturbs  the  bridge  negligibly. 
Calibrations  carried  out  for  helium  and  nitrogen  fit  the 
dependences  predicted  for  the  bridge  working  to  the  infinite 
impedance  within  2%  at  all  points.  Drift  of  the  system  is 
5  4%  of  full  scale  in  about  12  hours  for  liquid  helium 

sensitivity  and  about  ten  times  less  for  liquid  nitrogen 
sensitivitv.  The  stability  could  be  improved  by  the  use  of 
more  stable  operational  amplifiers. 


The  amplification  in  the  position  1  is  ten  times  higher  than 
in  the  position  2. 
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APPENDIX  VII 

DESCRIPTION  OF  A  TYPICAL  RUN 

After  the  glass  substrate  with  gold  contacts  is 
clamped  to  the  bottom  of  the  liquid  helium  pot,  the  copper 
shield  with  the  movable  mask  is  screwed  to  the  liquid 
helium  (  LHe  )  vessel.  The  mask  is  connected  by  piano  wire 
to  an  outside  dial  and  the  position  of  the  mask  is 
calibrated  in  terms  of  the  dial  reading.  Then  the  copper 
shield  is  fastened  to  the  liquid  nitrogen  (  LN  )  vessel  and 
a  movable  shutter  is  also  connected  by  steel  piano  wire  to 
the  outside  knob.  The  evaporation  sources  are  filled  with 

fresh  metals  and  the  system  is  sealed  off  and  pumped  (  for 

-7 

about  15  hours  )  to  a  vacuum  of  about  2  x  10  Torr.  Then 
all  evaporation  sources  are  covered  by  the  shutter  and  are 
heated  up  until  the  evaporation  starts,  thus  melting  and 
degassing  the  materials  which  are  to  be  evaporated  later. 
The  LN  is  passed  through  the  Meissner  coil  inside  the 

evaporation  space.  The  preparation  of  the  junctions  is 

-7 

started  when  the  pressure  drops  below  1  x  10  Torr. 

If  a  SiO  layer  is  wanted  on  the  substrate,  this  is 
evaporated  first,  then  the  mask  moved  and  the  aluminum 
base  layer  is  evaporated.  The  aluminum  thickness  is 
monitored  by  the  quartz  crystal  monitor  and  the  resistivity 
of  the  evaporated  film  is  measured  simultaneously.  When 


. 
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the  evaporation  of  the  base  layer  is  completed,  the  main 
gate  valve  is  shut  off,  the  flow  of  LN  through  the  Meissner 
coil  is  interrupted  and  oxygen  at  a  pressure  of  about 
0.1  Torr  is  introduced  into  the  system.  The  oxygen  is 
purified  by  passing  it  through  a  coil  immersed  in  LN.  Then 
the  whole  strip  of  aluminum  is  oxidized  using  the  glow 

discharge.  After  that,  the  system  is  pumped  to  a  pressure  of 

-7 

~  5  x  10  Torr  within  a  few  minutes.  Then  the  flow  of  LN 
through  the  Meissner  coil  is  again  initiated,  thus 

_  n 

improving  the  vacuum  pressure  to  ~  1  x  10  Torr.  The  mask 
is  moved  to  the  next  position,  pressed  against  the  glass 
substrate  and  the  first  layer  of  the  proximity  sandwich  is 
evaporated.  The  thin  layer  of  oxide  is  allowed  to  grow 
(or  a  layer  of  SiO  less  than  20  R  thick  is  evaporated)  and 
the  junction  is  completed  when  the  second  layer  of  the 
proximity  sandwich  is  evaporated.  The  same  procedure  is 
repeated  for  the  other  junction.  When  this  is  done  the 
shutters  on  both  LN  and  LHe  shields  are  closed.  The  flow 

of  LN  through  the  Meissner  coil  is  then  stopped  and  the 

-5 

vacuum  slowly  deteriorates  to  a  value  of  -  3  x  10  Torr 
within  10  minutes.  After  that  the  pressure  quickly  drops 

_7 

to  the  10  Torr  range.  This  effect  is  caused  by  gas 
molecules  trapped  on  the  cold  surface  of  the  Meissner  coil, 
which  are  released  as  the  coil  warms. 

After  the  junctions  are  prepared  the  level  detectors 

are  switched  on  and  all  three  vessels  are  filled  with  LN. 
(Nitrogen  vessel , helium  vessel  and  helium  pot, into  which  the  LN 


- 


« 
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is  syphoned  through  the  1/8  inch  diameter  pumping  hole  ) . 

The  filling  of  the  system  with  LN  can  be  finished  within  half 
an  hour.  Thermal  equilibrium  with  LN  in  the  nitrogen  vessel 
(  no  violent  boiling)  is  usually  reached  after  1  hour;  in 
the  other  two  vessels  the  equilibricum  is  reached  much  sooner. 
At  this  stage  it  is  safe  to  leave  the  system  for  about 
20  hours  without  refilling. 

In  order  to  cool  the  sample  to  LHe  temperature,  we 

proceed  as  follows.  A  thin  stainless  steel  tube  is  inserted 

into  the  helium  pot  and  LN  is  blown  out.  Then  the  helium 

-2 

pot  is  pumped  to  about  10  Torr  and  LN  is  blown  out  of  the  LHe 
vessel.  The  helium  pot  (  and  all  system  with  which  it  is 
interconnected  )  is  filled  with  He  gas  at  atmospheric 
pressure.  A  copper  rod  is  inserted  through  the  1/8  inch 
hole  into  the  helium  pot  to  provide  heat  contact  with  the 
LHe  vessel  (  through  the  exchange  gas  in  the  helium  pot 
system  ).  Then  LHe  is  transferred  into  the  LHe  vessel. 

After  the  temperature  of  the  helium  pot  reaches  about 
5  to  8°K  (  in  about  20  or  30  minutes  ) ,  the  copper  rod 
is  replaced  by  a  thin  syphon  and  the  helium  pot  is  filled 
with  LHe.  The  helium  level  is  detected  by  a  carbon 
resistor  placed  in  the  pumping  line  about  7  inches  above 
the  top  of  the  helium  pot.  At  the  beginning  of  the  LHe 
jgns f er  the  temperature  of  the  helium  pot  is  increased 
to  about  10 °K  by  warm  gas  from  the  inside  of  the  syphon 
tube.  About  10  litres  of  LHe  are  required  to  cool  the  system 
from  LN  temperature  and  fill  it  with  about  5  litres  of  LHe. 
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The  temperature  during  the  measurement  is  regulated 
by  controlling  the  vapor  pressure  of  LHe  in  the  helium  pot. 
The  lowest  attainable  temperature  is  1.05°K.  We  tried  to 
lower  the  temperature  by  pumping  both  the  helium  pot  and 
the  LHe  vessel  simultaneously;  however  the  temperature  was 
decreased  about  20m°K  only  at  the  lowest  attainable 
temperature . 

After  all  the  measurements  are  finished  the  LHe  is 
boiled  off  (  from  the  pumping  pot  by  a  heater  in  situ, 
from  the  LHe  vessel  by  a  dipstick  heater  )  and  the  system 
is  left  overnight  to  warm  up  by  radiation.  Next  morning 
the  main  gate  valve  is  shut  off  and  the  evaporator  - 
cryostat  is  filled  with  helium  gas  at  a  pressure  of  a  few 
Torr.  About  three  hours  later  the  system  is  all  at  room 
temperature . 
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